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Abstract: 

With a view towards implementation in microscopic transport simulations of heavy- 
ion collisions, the properties of spin-isospin modes are studied in nuclear matter con- 
sisting of nucleons and A isobars that interact by the exchange of vr and p mesons. 
For a standard p-wave interaction and an effective g' short-range interaction, the 
dispersion relations for the spin-isospin modes, and the associated amplitudes, are 
calculated at various nuclear densities and temperatures, within the random-phase 
approximation. Quantities of physical interest are then extracted, including the total 
and partial A decay widths and the A cross sections in the nuclear medium. The 
self-consistent inclusion of the A width has a strong effect on the A cross sections 
at twice normal nuclear density, as compared with the result of ignoring the width. 
Generally, the obtained quantities exhibit a strong density dependence, but are fairly 
insensitive to the temperature, at least up to T = 25 MeV. Finally, it is described 
how these in-medium effects may be consistently included into microscopic transport 
simulations of nuclear collisions, and the improvements over previous approaches are 
discussed. 
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1 Introduction 



In recent years, large efforts have been devoted to the understanding of hadronic 
matter at high densities and temperature, created in colhsions between two heavy 
nuclei at bombarding energies from a few hundred MeV per nucleon up to several GeV 
per nucleon. A large number of energetic particles are produced in such collisions and 
they may be used as probes of the hot and dense phase of the reaction, since they 
were not present initially 0, |^. 

Particle production in heavy-ion reactions has been fairly successfully described 
by microscopic transport models, such as BUU and QMD 0, |], |, In these trans- 
port models, the nucleons propagate in an effective one-body field while subject to 
direct two-body collisions. Sufficiently energetic nucleon-nucleon collisions may agi- 
tate one or both of the colliding nucleons to a nucleon resonance, especially A(1232), 
iV*(1440), and A^*(1535). Such resonances propagate in their own mean field and may 
collide with nucleons or other nucleon resonances as well. Furthermore, the nucleon 
resonances may decay by meson emission and these decay processes constitute the 
main mechanisms for the production of energetic mesons 

The by far most important nucleon resonance is the A(1232) isobar. It contributes 
substantially to the production of pions, kaons, anti-protons, and di-leptons, either 
directly via its decay, or indirectly as an intermediate particle [Q. The A isobar 
is also very important for the thermal equilibration of the reaction, by absorbing 
the kinetic energies of the nucleons. The transport descriptions normally employ the 
vacuum properties of the resonances and mesons, i.e. the needed cross sections, decay 
widths, and dispersion relations are taken according to their values in vacuum [Q. 
However, it is well known that the pion changes its dispersion relation in the nuclear 
medium, and that the A changes its decay width p. These in-medium properties 
have proven to be important to incorporate in the description of various reactions 
involving-isospin degrees of freedom, for example 7r-nucleus reactions [0, and (p,n) 



and ( He,t) charge-exchange reactions M, 10 



It is therefore reasonable to expect in-medium properties of vr mesons and A iso- 
bars to play an important role in transport descriptions of heavy-ion collisions. Some 
in-medium modifications have already been employed in calculations of heavy-ion 
collisions, both qualitatively [|T^] and by transport simulations |]T2|, [T^, . Although 
the particular properties considered in those studies may not change much in more 
refined treatments |T2|, ITB|, M , other properties can change significantly. More refined 



models have been used in various contexts, for example the self-consistent coupling 
of pions and A-hole states in nuclear matter [|l^ and the calculation of the A width 



in nuclear matter |T6|, . The purpose of the present investigation is to develop and 
employ such a refined model to derive several in-medium quantities that are useful 
for transport models, and to discuss their meaning and how they can be implemented 
consistently in dynamical simulations. 

It is important to recognize that the physical scenarios in heavy-ion collisions differ 
from those in charge-exchange and vr-nucleus reactions, and therefore the effects of in- 
medium properties may differ as well. Charge-exchange and vr-nucleus reactions occur 
in a cold medium at relatively low density (in the surface). As larger densities are 
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probed in heavy-ion collision, the in-medium effects are expected to be enhanced. On 
the other hand, there are other mechanisms that may reduce the in-medium effects: 
the dense medium is usually hot as well, many individual baryon-baryon collisions 
and resonance decays are involved, and many of the observed mesonic particles have 
been created at the surface of the reaction zone and are thus less sensitive to the 
dense interior region. 

In infinite nuclear matter, a system of interacting tt mesons, nucleons, and A 
isobars will couple to form spin-isospin modes. Some of these modes are collective 
and correspond to pion-like states (quasipions), while other modes are non-collective 
in their character. While a free pion has only a pion component in its wave function, a 
collective pionic mode contains components from many AA^-^ and A^A^-^ states. The 
strength of the various components will vary with the momentum of the collective 
mode, and the energy will differ substantially from the energy of the unperturbed 
pion or AN~^ states. On the other hand, non-collective modes are dominated by a 
single component of a baryon-hole state (iViV~^ or AA^~^) and its energy is close to 
that of the corresponding unperturbed baryon-hole state. 

Our aim is to calculate the properties of the spin-isospin modes in infinite nuclear 
matter at various nuclear densities, pn, and temperatures, T. From these properties, 
we will also deduce a number of other physical quantities, such as decay widths and 
cross sections, which will depend on pat and T. In order to obtain a consistent de- 
scription, it is important to calculate simultaneously the dispersion relation, the A 
width, and the cross sections. As new decay channels for the A are present in the 
nuclear medium, the A width will differ from its vacuum value|16, T7|. However, some 



of these decay channels correspond in a transport description to A-nucleon collisions 
and should be excluded from the A width in a transport description. Instead, it is 
important to take into account the in-medium properties for calculating A-nucleon 
cross sections. It is also important to calculate the dispersion relations and the A 
width simultaneously, because these quantities are interdependent. We have done 
this with an iterative procedure. In a subsequent paper we will present results from 
microscopic transport simulations of heavy-ion collision with the presently calculated 
in-medium properties having been implemented by means of a local density approxi- 
mation, Pn = PN{,f) and T = T(r). 

The effects of including the in-medium properties of pions and A isobars in descrip- 
tions of heavy-ion collisions have been studied previously but only some particular 
aspect was considered, such as the A production cross section , the production of 



pionic modes |TT|, or the changed dispersion relations of the pionic modes |[T^, 
rather than treating all of the aspects consistently within the same model (dispersion 
relations, A width, and cross sections), as is done in the present paper. Furthermore, 
we have improved upon a number of approximations done in previous works. While 



the A width was excluded from a number of quantities in refs. [|12|, |13|, |I4[, we have 



included the A width in a self-consistent manner.^ Some of the previous treatments 



^We include the A width by an iterative procedure and refer to this as the self-consistent case. 
This nomenclature is not intended to suggest that the entire set of coupled equations for the in- 
teracting hadons are solved self-consistently. Such truly self-consistent calculations were recently 
carried out with the A and tt degrees of freedom included, but ignoring the NN^^ excitations [p^. 
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excluded the nucleon-hole channel and used a simplified expression for describing the 
continuum of A-hole states [l^, and the in-medium properties were calculated 
in cold nuclear matter in refs. |T^, |1^, |13|, |T^ . The avoidance of these approximations 
has important consequences for the dispersion relations, the A width, and the cross 
sections. This will be further discussed in sec. ^, where we also will emphasize the 
importance of treating collective and non-collective spin-isospin modes in a consistent 
way within the model and will compare to previous treatments. 

Although our present study constitutes a more consistent way of incorporating 
in-medium effects into a dynamical transport description than what has been done 
previously, our treatment does rely on certain approximations. Thus, we limit our 
considerations to systems consisting of interacting nucleons, A isobars, and vr and p 
mesons, and we assume that only relatively few A isobars and mesons are present. In 
addition to the increasing role played by higher resonances, there may also be other 
effects occurring at high densities that are not taken into account in our model, such 
as partial restoration of chiral symmetry. 

In sec. we present the model. While part of the formalism can be found in 
the literature {e.g. ref. 0), there are also important differences from the traditional 
formalism. Therefore, to make the presentation as clear as possible, we include some 
steps that can be found elsewhere. In sec. |], we motivate and discuss our choice 
of parameter values. The results are then presented in sec. § which also contains a 
discussion of the implementation of the in- medium quantities in transport simulations, 
and a comparison with previous works. Finally, our results are summarized in sec. ^. 



2 The model 

We will consider a system of interacting nucleons (A^), delta isobars (A), pi mesons 
(tt), and rho mesons (p). In order to investigate the matter properties of the inter- 
acting particles, we employ a cubic box with side length L; the calculated properties 
are not sensitive to the actual size, so we need not take the limit L — oo explicitly. 

The in-medium properties are obtained by using the Green's function technique, 
starting from non-interacting hadrons. The non-interacting Hamiltonian can be writ- 
ten 

Ho = Y ^kbih + Y + ^^pMpiPn ■ (1) 

k I n 

Here the index k = {pf.; Sk,ms^] tk,mt^) represents the baryon momentum, spin, and 
isospin. The spin and isospin quantum numbers, Sk and tk, take the values | and | 
for N and A, respectively. The energy of baryon k moving in a (spatially constant) 
potential is denoted e^. The baryon creation and annihilation operators, bl and 6^, 
are normalized such that they satisfy the usual anti-commutation relation, 

{bl h'} = ■ (2) 

In the pion part of Hq, the index / represents the pion momentum and isospin, / = 
{Pi,Xi = 0, ± 1), while the index n in the p term also includes the spin ms„. The 
meson energies are given by hu^^^p = [rri^ p + q^Y^"^ respectively, and the creation and 
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annihilation operators of the pion are normahzed such that they satisfy the usual 
commutation relation, 

[^/) ^v] = kv ' (3) 

and analogously for the p meson operators, p\ and 

Note that the A isobar described by Hq has no decay width. Fa = 0. When the 
interactions are turned on, the A width will emerge and it will then automatically 
include also the free width. 



2.1 Basic interactions 

At the NtcN and A^ttA vertices we will use effective p-wave interactions, Vnt^n and 



VnttA, which in the momentum representation can be written as 

1 

2mArC^ 



Vi 



NttN 



Vj 



NttA 



IC 



IC 



M 

{hey 

L3 



mAC^ + y/s 



NN 



2 f TT 

J NA 



FAq) i^-QcJ r • 0, (q) (4) 
F^q) -i^ (q)+h.c.(5) 



In these expressions, y/s is the center-of-mass energy in the Nn system and q^^ is 
the pion momentum in the A^vr center-of-mass system, which in the non-relativistic 
limit is given by 



rriNC 



q- 



Pn 



rriNC^ + huj ^ rriMC^ + huj ^'''^ ' ^ ^ 

where hu and q is the pion energy and momentum, and is the nucleon momentum 
in an arbitrary frame. The Pauli spin and isospin matrices are denoted cr and r, and 
and are spin and isospin ^ to | transition operators normalized such that 
•Si 1^, H >= 1-^ The momentum representation of the pion field is taken as 



< 



3 3| 

2' 2l'-^+ll2' 2 



^l{q) 



i2huAq) 

The interactions contain a monopole form factor. 



-q) 



(7) 



A? 



Al - (cqy 



(8) 



and the coupling constants are determined at (cg)^ = {huY — (cq)"^ = (m^c^)^ 
and ^/s = tjinc^ or y/s = rriAC^. Apart from the factor 2mBC^ /[mBC^ + a/s] the 



^For clarity, we generally employ bold- face characters to denote quantities with vector and tensor 
properties under ordinary spatial rotations, while arrows are employed to indicate the transformation 
properties under rotations in isospace. 
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interactions Vnttn and VatttA are obtained in the non-relativistic limit from the phe- 



nomenological Lagrangians [18 



m-,rC 



(9) 
(10) 



The factor 2mBC^ /[mBC^ + \/s] is a relativistic correction that takes into account 
that in relativistic calculations the energy denominators usually appear in the form 
2mBC^/[ (msc^)^ — s], while in non-relativistic calculations the form [m^c^ — y/s]~^ 
usually appears With this correction, the correct relativistic form of a Breit- 

Wigner resonance is obtained for the free A resonance [^. 

The interactions at p meson vertices are less well determined than for tt meson 
vertices. Here we choose a form analogous to eqs. (|) and (||), 



Vr 



NpN 



NpA 



IC 



{hey. 

{hey 
~1? 



NN 



nipC 



^F,{q) (o-xq,J.[?.</>,(g)] (11^ 



nipC 



+ h.c. (12) 
Apart from the factor 2mBC^ /[mB(? + these interactions can also be obtained 



as the non-relativistic limit of relativistic Lagrangians |l8 



In addition we will also include effective short-range interactions at nucleon-hole 
vertices, again written in momentum space, 

3 / r-^ \2 

,2._ (13) 



V, 



NN-^,NN- 



9nn 



In 



NN 



\Fg{q)\' {CT, ■ ■ T2) 



and the corresponding interactions obtained when one (or two) of the nucleons is 
replaced by a A. The strength of the short-range interactions is determined by the 
correlation parameters g'j-^j^f, q'mai S'aa- form factor 



A2 

9 



Al - {cqf 



(14) 



is omitted, this interaction has vanishing spatial range, V ~ 5{ri — 



2.2 Green's functions 

The interactions between A^, A, vr and p will lead to the formation of spin-isospin 
modes which will carry the quantum numbers of the vr or p mesons. The propagation 
of these spin-isospin modes can be represented by an appropriate Green's function. 
For this purpose, we define a Green's function 

iG{a,p-t,t') = l< Q{Bl{t)BJt')} > , (15) 
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where G is a time ordering operator, and the brackets < ■ > denote either the 
expectation value of the interacting ground state (at zero temperature) or the thermal 
average (at finite temperatures). The interacting ground state is an eigenstate of the 
Hamiltonian, if, which consists of the noninteracting part Hq (given in eq. (|l|)) and 
an interacting part defined by the interactions in eqs. (^|-[T3|). The "channel" indices 
a and (3 are used for convenience, so that the creation operator 

represents either a vr meson operator, 

^i = ^L(9a), or Tf_^^{-qJ (17) 
a p meson operator (defined analogously), or a two-baryon operator, 

Bi = > (18) 

which in the zero-temperature limit corresponds to creation of a NN~'^ or AA^~^ 
state. 



At zero temperature the Green's function in eq. ( |T5|) can be written in the energy 
representation as 



iG{a,P;iu) = / dit-t') e''^^'-''UGia,P;t,t') 

J —oo 



huj — TiUy + irj 
huj + huJi, — irj 



^ hid — hujy + irj ^ huj + huJiy — irj 

where, as is usual, we have approximated the excited states {"^i, > by excited RPA 
states, > ~ Ql l^^o >, generated by a generalized RPA operator, 

Ql = j:x:Bi, (20) 

a 

and approximated the matrix elements by 

< ^o\Bi\^. > = < ^o\[BlQl]\^o >^< %\[BlQl]\% > , (21) 
with |$o > denoting the non-interacting ground state. 



2.3 RPA approximation 

We want to calculate the spin-isospin mode Green's function defined in equation (|T3|) 
within the RPA approximation, symbolically 

G^P^(a, P; Lu) = G'.ia, /3; ^) + E G[,(a, 7; u;) V(7, lu) G^p^(«:, (3; u) . (22) 

7,K 
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The spin-isospin modes, here represented by the Green's function G^^^, will in this 
approximation be obtained as an infinite iteration of (non-interacting) pion, p meson, 
nucleon-hole, and A-hole states, represented by the diagonal Green's function G'q, 
coupled with the interactions specified in eqs. (|^-[13D which here are summarized by 
the symbolical interaction V. G^^^ is graphically illustrated in fig. |1|. As the lowest- 
order Green's function we will, in the case a corresponds to a meson index, take 



G'^ = GY or G'^ = G^'^ with 



G^^P{a,f3;u}) 



1 + 



+ 



n. 



7T,p 



huj — huJTr,p + iv Jiuj + huj. 



IT] 



n. 



n. 



7r,p 



huj + hcUj^n — if] huo — TiUJ-j^n — if] 



Sal3 


(23) 




(24) 




(25) 



Note that the free tt or p meson propagator, D^, is related to G^'^ and Gl^^ by 

1 



[GY{o^,P) + G'!:'{a,P)] . 



(26) 



In nuclear collisions at beam energies up to about one GeV per nucleon, which 
is the domain of application that we have in mind, only relatively few mesons and 
isobars are produced and so the associated quantum-statistical effects may be ignored. 
Accordingly, we assume ~ 0, ~ 0, and np ^ 0. Since we consider thermal 
equilibrium at a specified temperature T, the nucleon occupation probabilities are 



nN{k) 



1 + e{ek~f^)/T 



(27) 



In the case a corresponds to a two-baryon index we take G'q = G'^^ " or Gq 



G 



AN- 



with 



G 



NN- 



G 



AN- 



nN{k ) - UNih 



Sat 



hu - ei^ + Ck^ + ir] ■ sign(ci;) 

nN{ka) 

hu - ei^ + Cfc^ - SAAr-i(efc^ + LJ, Q + ir] 

riNila) 



(28) 



6n 



hu + Cfc^ - e/^ + Saat-i (e/^ - cj, fc^) - ir] "'^ ' 
The Green's function G^^ ^ has been calculated from 



G 



AN- 



dE 



(29) 



(30) 



Sal3 - I — 



27r 



G'^iL^E + u) G^ (ka, E) + G'^ {la, E + uj) G^{ka,E) 



N I 



where G^ is the full in-medium Green's function for the A, containing the A self 
energy Sa- We note that the quantity Saat-i in eq. ( pQ]) is identical to the A self- 
energy Sa when G^^ is calculated from eq. (pTI). However, as a first approximation 
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when calculating G^^^, the quantity Saat-i in eq. ( p9D can be ignored. The G^^^ 
obtained with this approximation can then in turn be used to calculate the A self- 
energy. It is therefore convenient to use different notations in order to distinguish 
between the input T,an-^ in eq. (|^) and the calculated A self-energy Ea using G^^^. 

We will obtain two different, but equivalent, expressions for the Green's function 
G^^^, as defined by eq. (|22D . The first expression is based on a summation of the 
Green's functions in eqs. according to the diagrams in fig. |I[ The details of 

this derivation can be found in the literature, for example in refs. 0, so in this 
paper we will only state the final expression in sec. p2.4| . This expression is useful for 
calculating quantities hke the total A width and different cross sections, involving a 



A isobar, and will be used in section 



The second expression is based on an expansion in RPA eigenstates. For this 
purpose we will derive a set of RPA equations, equivalent to eq. (^21) . We will solve 
these equations to obtain eigenvectors and eigenenergies for the different spin-isospin 
modes. The eigenvectors will yield the amplitudes of the different components (vr, p, 
NN~^, AN^^) forming the particular spin-isospin mode with the given eigenenergy. 
We will show that the eigenvectors form a complete orthogonal set, and we wiU 
expand the Green's function G^^^{a, P;u) on this set. This expression for G^^^ will 
be useful for calculating partial contributions to the total A width. Furthermore 
the RPA amplitudes of the different components will contain important information 
about the nature of the different spin-isospin modes. This will be further discussed 



in sec. 4.1 



2.3.1 Interactions and operators 

It is convenient to rewrite the total Hamiltonian specified by eqs. (P^-pT^) in the form 

H = H'o + + v^(4) ^ (31) 



+\ E Wb'M' ' (32) 

jkj'k' 

^^'^ = E^k' ^]^A' + E^.V ^]Pfc^' + h-c. , (33) 

jkj' jkj' 

= ^ E ^j>^fk' Wkk'h' 

jkj'k' 

= ^ E ""jkj'k' b]bkbk'bj' + i E ^jkj'k' b]bkbk'bj' ■ (34) 

jkj'k' jkj'k' 

The interaction V^^^ corresponds to the baryon-meson-baryon vertices in eqs. (^), 
(^, (|TT|) , and (|12D. In the spin-isospin summation the spin- longitudinal and the spin- 
transverse channels are orthogonal and do not mix. The interaction V^^^ contains 
the effective short-range g' interaction ([131) and is separated into a spin-longitudinal 
part and a spin-transverse part. In addition to these interactions, we also include an 
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effective interaction, v^J^,/^, , wliicli for tlie AA^~^ states includes tlie A self energy in 
our RPA formalism in a way analogous to eqs. (|22|) and (^). In the above expressions, 
Vjkj'k' = Vjkj'k' — Vjkk'j' is formally an anti-symmetrized two-body matrix element, but 
the exchange term vjkk'j' will be neglected in the calculations. 

The spin-isospin modes (or excited RPA states) |\E'^ > are created by an operator 
Ql- For specified values of their momentum q, isospin A, and spin fi the associated 
energy can be obtained as a solution to the RPA equations. In the spin-longitudinal 
channel the spin is zero, while the spin-transverse channel has spin 1, with two non- 
vanishing contributions, fi = ±1, for the spin-projection along the q-axis. We take 



q,A) 



(35) 



jk 



and 



Qliq, A, /i) = 5: X*,(g, A, fx)b]k + E Zliq, A, /i) Pk^ - E W^q, A, p) p, . (36) 

jk k k 

In appendix |^ we will restrict the summation over baryon and meson states in eqs. 
(H) and (H), by taking Xjk oc 5p^^p^+q, Zk oc 5p^,q5x^,x, and Wk oc 5p^_qSx^:-x. 
The RPA equations are obtained from the relation 



< [6Q, [H, gt]] >= hu < [6Q, > 



(37) 



with 5Q = blpj, TTj., Til, pr, or pj, and where the brackets < ■ > as previously denote 
the thermal average (which at temperature zero becomes the expectation value in the 
interacting ground state). 

2.3.2 RPA equations 

The general structure of the RPA equations is similar in the spin-longitudinal and 
spin-transverse channels. We will therefore not write out the symbols / and t in this 
section. Calculating the necessary commutation relations using < >= n{k), we 
obtain 



/ ^(1) + A^2) j^^ 

Ct D 

V -c^ 





(3^ 



with 



^(1) 

jkj'k' 

.(2) 
jkj'k' 

C jkr 

D„.' 
Afjkj'k' 



[cji — Ck' + T,jik']5jij5kik 



Lt 



[n{k) - n{j)]Vjk,i,j, 

Vjrk 

huj.j^(^qr^6rr' 

[n{k') - n{j')]Sj>jSk'k 



(39) 

(40) 
(41) 
(42) 
(43) 
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1 t- 

2' ''J 

1 t- 

2 ' J 

^ t- 

2' J 



(44) 



where the matrix C has the properties Cjkr = Ckjr and Cjkr = — Cj^r after the 
spin-isospin summations have been performed. Because we have 

[0 ih-' ' - ' 

Ejfc = < Eaat-i (cfe + hu, Pf. + q) {tk 

[ -SA7V-i(ej - - q) (4 - 2' -^j - 2 

the matrix A^^^ is non-Hermitian and so the usual RPA orthonormahty relations do 
not hold. In order to construct an orthonormal set, we will use the solution of the 
equation obtained by replacing A*^^^ by (A*^^))* . This corresponds to employing a bi- 
orthonormal set. The solutions of the auxiliary equation, as well as other quantities 
related to the auxiliary equation, will be denoted by a tilde symbol. 



, W ] . (45) 

58f) then 



The RPA equations have the property that if uoj^ is a solution of eq. 
uju = is a solution of the auxiliary equation. The equations (|38|) also have the 
property that if {NX" , , W")'^ is a solution of (|38D with the eigenvalue ujy^ then 



[UX%i{lo^) 



/ \MX-\k{uJ, 
V Z'^iio,) 



(46) 



is a solution of ( ^81) with the eigenvalue tu^ = —u^. 

Taking into account the strong cu dependence of S^at-i in the matrix A^^\ but 
neglecting other weak u dependences due to form factors and relativistic corrections, 
we obtain the orthonormahty relation 



6,^ sign(Re u,) = ( X^^, Z\ 




(47) 

Here the u dependence of S^^-i is contained in the factor Vjkjik'itjtk,tj'tk'), which 
has the following non- vanishing elements 



--) 

'ijk,jk\2 2 ' 2 2 
,31 31. 
•22' 22' 



Vjk,jk\2 2 ' 2 2 



1 - 



SAAf-i(e,- 



(48) 
(49) 



SAAf-i(efc + huJf^pA - SAAr-i(efc + huj^,pj 



htu* 



(50) 



While the u dependence of Faat-i = —21m Saat-i is strong, the u dependence of Re 
Eaat-i is weak and can be neglected. It is convenient to take the quantity r]^''^ in the 
approximate form r]jj!'jj^ ^ [VjkYVjk 'with 



[1 + idT/^N-^ (e/c + h^^u, Pk + q)/'2dhuj 
1 + idT/^]^-i{ej — hujy,p, — q)/2dhuj 



1/2 
1/2 



tk 
tk 

tk 



2' 3 
2' ''J 



(51) 
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and fij,. = {rjj^)*- With this approximation, the resolution of unity takes the form 



sign (Re 



( r]^X' 



-[W 



and an arbitrary vector F can be expanded as 



(52) 




(53) 



We have checked numerically that the orthonormality and completeness relations, 
( ^71) and (|52D , with the approximation (|5TD , are satisfied within an error less than 
3-4%. 

It is convenient to express the two-body Green's function as 





(M 







(M 







G^P^(a;) = 





1 


j Giuj) 





1 






^ 







^ 








(54) 



where G{uj) can be expanded on the RPA states as 



(55) 



When T]'^ is unityi we can use eq. (|^) to determine the solution, 

sign (Re u^) 



huj — hojy + i5 sign(Re u^,) 



(56) 



After rewriting the expansion of G^PA as a sum over positive frequencies by using eq. 
(^6]), we arrive at the following expression for G^^"^, 



Re u)y>0 



huO — hujy + 17] hu + hUy — 17] 

where ?/> and ?/< are short notations for 



(57) 



and 



[AfX'^M]i^,,^ 

z^M 



{51 



At T = this is the same expansion as in eq. (|T^), with the approximation in eq. 
(pTI). However, the expansion in eq. (|57D also holds at finite temperatures. 



•^This is generally the case when Faat-i is independent of lo, and in particular when Faat- 
vanishes identically, as in our reference case. 
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However, when rj'^ is taken from eq. (51), with dT/^j^~i/dhuj ^ 0, the expansion 



of G^^^ will no longer be diagonal in the channels z/, /z, i.e. g'-''^{uj) ^ g^{uj)5y^. 
Therefore, in the present paper, the RPA expansion of G^^^ will only be used when 
Faat-i vanishes. 

It is straightforward to obtain the explicit solution of the RPA equations for 
the interaction specified by eqs. (^-[131), but the notation is somewhat tedious. We 
therefore present the essential steps together with the final expressions in appendix 

2.4 The A self energy and A cross sections 

In this section we will derive expressions for the A self energy Sa and the cross 
sections (t(1 + 2 — 3 + 4) where 1-4 represent baryons. These quantities can be 
expressed by an effective interaction that is obtained by inserting G^^^ between two- 
baryon states, as illustrated in fig. 0. We will denote this effective interaction with 
M(34, 12). It is convenient to write the spin-isospin matrix elements as a separate 
factor and define the quantity M''*(34, 12) by 

M''*(34, 12) = ^9''*(31)(^9''*(24))*M''*(34, 12) , (59) 

where the spin-isospin matrix elements, with l=2=4=iV and 3= A, in the spin- 
longitudinal (/) and the spin-transverse [t] channels are written 

^\'il){'d\2A)y = <m,3|(5+-qr3i)|m,, ><m,J(o-.qr24)l"^.2 > 

<mts\T \mt^ > • < mt J f Im^^ > (60) 

#(31)(^*(24))* = <m,3|(5+xg3i)|m,, >.<m,J(o-xq24)l^.2> 

<mt.^\T \mt^ > • < mt^\ T \mt^ > . (61) 

When G^^^ is expressed as a sum of non-interacting Green's functions, the effec- 
tive spin-isospin interaction M(34, 12) becomes 

M'(34,12) = -p^AD^Flql^{UA2) + F^^g'^^{U,l2)\ (62) 

M*(34,12) = (34, 12) + -^i^^ (63) 

The expression of the quantities -D7r,p, g^ff(34, 12) and (7gg(34, 12) are somewhat lengthy 
and are therefore given in appendix |B[ 

Alternatively M(34, 12) can also be expressed using the RPA expansion in eq. 

\ /i''*(31; v)h^^\2A] v) /i''*(31; z/)/i''*(24; z/) 1 



Re aj">0 



M^'*(34,12) = 

l.t 

Jl 



huo — huhf + it] hu + huii\f — it] 



fvr -fn 

+ TZrZ^^Q 934,12 ■ (64) 
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The factor h{jk, v) is obtained from the interactions at the vertex consisting of 
baryons j and k, and the spin-isospin mode v. The interactions to be used de- 
pend on the non-interacting states that the mode is made up of and must therefore 
be multiphed by the amphtude of the corresponding state. For example, 

h\3K i^)^\3k) = [Z\u) + W\u)] + E V,u,mn , (65) 

V ^'it^TT mn 

where V^vfc is defined in (^) and (|^), Vjk^mn is defined in (]TB|), and the amphtudes 
X^^, Z\ are defined in eq. (^). From appendix ^ we obtain the amphtudes 
{Xmn, Z, W) as the solution of the RPA equations. These expressions are also some- 
what lengthy and they are therefore relegated to appendix 0. 



2.4.1 The total A width 

The A self energy S^, is calculated according to the diagrams in fig. ^, by taking into 
account all the diagrams corresponding to the A decaying into a spin-isospin mode 
and a nucleon, which then again form a A. Since the spin-longitudinal and spin- 
transverse channels are orthogonal we can treat them separately. The spin-isospin 
summation gives, for a spin average over the external A spin states, a factor 1/3 in 
the spin-longitudinal channel, and a factor 2/3 in the transverse channel. 

The contribution to the A self energy in the spin-longitudinal channel can be 
expressed as (see also ref. [|1^, 0), 

3 



Sk(i?A,PA) = M T I E r ^ M'(AiV,iVA) G^(Ea - nu^p^ - q) . (66) 



By writing the nucleon propagator Gat as 



E-eN{pN)+ir] E - eiv(p^) - zr/ 

we can carry out the uj integration in (|66|) by performing a Wick-rotation |[T6| , |19|, 
and we obtain the expressions for the A width. Fa = — 2Im Ea, as 

F^(EA,PA) = Im^(^) Y.m£)-n{p^-q)]M\AN,NA) . (68) 

where the energy available for the spin-isospin mode is given by 

S = E^-cn^Pa- q) ■ (69) 
In the same way we obtain the contribution to Fa from the transverse channel. 



4 /fir\ ^ 

F^(Ea, Pa) = Im - l^-j E " ^(Pa " q)] M\AN, NA) 



(70) 



13 



2.4.2 Specific A channels 

The total A width, Ta = + T*^, gives the transition probability per unit time for 
the A resonance to decay to any of its decay channels. In a transport description one 
explicitly allows the A resonance to decay into specific final particles. Consequently, 
one needs not only the total A width (which is the sum of all decay channels) but 
also the partial widths governing the decay into specific RPA channels. These decay 
channels consists of a nucleon and one of the spin-isospin modes. Since we have 
access to all the amplitudes of a given spin-isopsin mode on the different unperturbed 
states, it is possible to derive an expression for the partial contribution to Fa from 
the A decay to a specific mode v. The right-hand side of fig. ^ shows a diagrammatic 
representation of such a process. The partial A width for a A decay to a nucleon and 
a spin-longitudinal mode v becomes 

x^iv(PAr) (27r)^5(pA -Pn- q)'2n5{EA -Cn - hu) 
= MrJ -^,\h\AN,uj,)\^nN{pA-q)27TS{EA-eN-nuj,) (71) 

where un = 1 — it-n, and the spin-isospin summation gives a factor A4r = 1/3 in the 
spin longitudinal channel and 2/3 in the spin-transverse. This expression is identical 
to the contribution from one of the u terms in eq. (|68D, if the RPA form ( [6^ ) is used 
for M(34, 12) and Im uj^, = 0. If the Faat-i is taken to be zero in the calculation 
of the spin-isospin modes the energies Ui, will be real. In this case the summation 
over modes u in eq. (|64D corresponds to a summation over physical decay modes, 
A ^ N + u, when the A width is calculated from (|68D . 

When the A width is included self-consistently in the calculation of the spin-isospin 
modes the energies will be complex, Im ct;,^ < 0. This implies that the energy of 
the spin-isospin mode u no longer is distinct, but instead has a Breit-Wigner-like 
distribution with a width 2 Im hui, centered around Re hu^. To obtain the partial 
decay width to the mode u, F^, we therefore need to "sum" over all possible energies 
of the mode u, using that the probability to find the mode v in the energy range 
between e and e + is given by the factor 

' ^"^^^^ de. (72) 



TT (e — Re TiujyY + (I^i huj^ 
The expression for the partial A width is thus modified to 

F^(i?A,PA) = Mr j -^,de\h\AN,e)\^l-nM{PA-q)] 

1 Im hu^ n x/n^ \ 

^ ~1 5"^^ — y^VT^ — I — 27r(5 Ea - e^v - e 

vr (e — Re nuOyY + (im TiuOy)'^ 

d^q 
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X 



Im hui, 
(e — Re huJuY + ^^u) 



(73) 



Unfortunately the expression ( [75D cannot be used directly since we do not calculate 
the amplitudes of different modes with real energy e, but rather the amplitudes of 
a single mode having a complex energy hu,^. Therefore the amplitudes X^^, 
and W" will be complex quantities and thus also h{AN,u,^). However, the squared 
amplitudes 

VmnXUVmnKJ* , Z^iZr , and W^iWr 

have very small imaginary parts. We can therefore obtain a good approximation for 
h{AN, e) by taking 



[C\u) + w\u)]+J2V,k,mnCni'^) , (74) 



with 



c 



Re n'^ 'X'^ (ri'^ \ 

Imn rnn\ Imn i 
, ~ X 1 1/2 

Re ^'^(ly^^*^ 



mnl 



1/2 



Re Wiy^' 



1/2 



(75) 
(76) 
(77) 



We have numerically compared the total width obtained by summing the partial 
widths based on this approximation, with the correct total width based on eqs. (|6^), 
(H), (HD, and (|70D. We find, in the range of invariant A masses m from 1000 MeV/c^ 
to 1400 MeV/c^, that the approximation of the partial widths somewhat over-predicts 
the total width at low invariant masses and somewhat under-predicts it at large 
invariant mass. The relative error is between 0% and 20% depending on m. To 
improve our approximation we therefore multiply the approximate partial widths! by 
a factor ca(-Ea;Pa)5 depending only on the A energy and momentum, to obtain the 
correct total width when all partial widths are summed over. The uncertainty of the 
partial widths obtained by this procedure should thus be quite small. 

In a transport description where collective modes are propagated as quasiparticles 
one will need a cross section for the inverse process z/ + — > A. This cross section 
is obtained analogously to the partial A width in (^). When a A with mass m is 
created in the z/ + collision the cross section is written 



(T(z/ + Ar^ A) = [Ticf M,N\h^'\^X,^,)\ 



2 m(? m^(? 



47rp3(m2 



(7J 



where the factor p^{m?), defined in eq. (|83D , takes into account the finite width of the 
A, and the spin-isospin factor M-^n is 2/3 in the spin-longitudinal channel and 4/3 
in the spin-transverse channel for the process n^'+p^ A++. 



"^The partial width is not multiplied by ca since the amplitudes and energies of the 

nucleon-hole modes are real and do not depend on Faat-i- 
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2.4.3 Inelastic nucleon-nucleon and A-nucleon cross sections 



Next we will derive the formulas for calculating cross sections for the processes 



N + N ^ A + N and A^ + A^A^ + A^. 



(79) 



We will start by writing down the 5- matrix for the processes of interest. The pro- 
cedure is very similar to the procedure for deriving Sa- We will first present the 
formalism in the spin-longitudinal channel and then generalize to the spin-transverse 
channel. 

In what we will refer to as the direct term, we consider baryon 1 colliding with 
baryon 2. Baryon 1 will after the collision appear as baryon 3, while the incoming 
baryon 2 becomes 4 after the collision, see fig. ^ Baryon 2 may be either a or a 
A depending on the process, and in the same way may baryon 3 be a A or a A^. We 
denote the transferred energy and momentum hj ud and q^, respectively, with 



hut 



D 



63 - ei = 62 - 64, 



Qd=P3-Pi=P2- Pa 



^0) 



We also take into account the exchange process where baryons 3 and 4 are inter- 
changed. 

The cross section is obtained from 



where the transition probability per unit time is given by 



^1^ 



1^ 



w 



At 



m 



Sfi is the total scattering matrix of the process and At is here a finite time interval 
which will tend to infinity at the end of the calculation. The function p^{m1) takes 
into account the finite width of the A when baryon 3 is a A, 



m 



TT (m^-m^)2+m^rA("i.3)^ 



3 
3 



N 
A 



^3) 



In contrast to the calculation of the A self-energy the spin-isospin summation is 
not performed in the S-matrix, so the spin-isospin matrix elements need to be kept. 
In addition, an overall energy and momentum conserving (5-function is included, and a 
factor ^rriiC^/ei for each external baryon. We also get a contribution from the short- 
range interaction, eq. (|13D, acting directly between the vertices 3 <— 1 and 4^2. 



Taking this into account we can write down the contribution from the direct term in 
the spin longitudinal channel to the S'-matrix 



s^n = M\^ (34,12) n 

V -'^ / 7=1,4 L . 



2Tc6{ei + 62 - eg - 64) 5p^+p.-,,p,^+p^ ■ (84) 
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The contribution from the exchange term is obtained by interchanging baryon 3 and 
4, which gives 



hujE = e4-ei = e2-e3, Qe = Pa - Pi = P2 - Ps (85) 

The expression for is identical to S'^i with the replacements: lud — > ue, Qd ~^ Qe^ 
31 ^ 41 and 24 ^ 23. 

The total ^-matrix can now be written 



\ 3 - ^ 

he 



with 



rrijC 



27r(5(ei + 62 - eg - 64) Sp^+p^ p^+p^ , (86) 



M = M'(34, 12) + M'(43, 12) + M*(34, 12) + M*(43, 12) , (87) 

where the expressions for M* in the spin-transverse channel are obtained analogously. 

Summing over final states and averaging over initial states, we obtain the differ- 
ential cross section in the center-of-mass system as 

1^ - (*''^2^2sVt^ /*"^ S |M|^( n ,„,c>3(ml) . (88) 

spm J-i'"' 

From this differential cross section we can also write down an invariant cross section, 
= {hcf^ ^— [ dml y \M\\ W m.^)p3{ml) . (89) 



spm 



In these expressions we have used that the relative velocity Vrei can be be expressed 
by means of the relativistic invariant 

/ = exe^Vreilc = \l{eie2 - cp^ ■ cp^Y - {mic^Y{m2C^Y . (90) 



3 Parameter values 

The A-hole model with vr and p meson exchange and an effective g' short-range in- 
teraction contains a number of parameters. Although these can presently not be 
determined uniquely, it is possible to use existing experimental information to put 
a number of constraints on the parameters, thus limiting the range of their values. 
In this section we will present our choice of parameter values together with a discus- 
sion and a motivation of this choice. For convenience, all the parameter values are 
summarized in table |l|. 

For the tt and p coupling constants, f^^^ Ina^ Inn^ /at a 5 choose their values 
such that they are consistent with vr-absorption data on the deuteron, i. e. according to 
the range of possible values in ref . |2^ . For the tt coupling constants this together with 



the condition that our model reproduces the value of the free A width at resonance, 
r|^^(mA = 1232 MeV/c^) ^ 115 - 120 MeV , (91) 



17 



fixes the value of f^j^ = 1.0 and f^^ = 2.2. For the p couphng constants we 
follow ref. and relate /j^^ and /j^^ according to the quark model relation, f^^ = 



72/25/^^. However, the value of is not fixed by the work in ref. pOl, but is only 
constrained to be in an interval of possible values, approximately 5.6 < < 7.8. 
In our work we choose the value /j^^ = 6.2 from this interval. This value has also 
been used in other studies related to the present work, e.g. in ref. [|12|. With this 
choice we are able to reproduce free cross section for the process p + p ^ + n, 
if we adjust the remaining parameters (essentially g') appropriately. This is not the 
only choice of that can reproduce these cross sections, but for other choices of 
the remaining parameters would have to be readjusted. 
Closely related to the values of the coupling constants are the values of the cut-off 
factors, and Ap, in the monopole form factors that are included in the interactions 



used at the vertices with a vr or p meson. In ref. pO|] the values A^, = 1.2 GeV and 
Ap > 1.5 GeV were used. In our work we take a somewhat lower value of the pion 
cut off factor, A^r = 1.0 GeV, and we take Ap = 1.5 GeV, to achieve a relatively fast 
cut off in the numerical integrations. 

When we calculate the dispersion relations, i.e. find the energy-momentum rela- 
tion huj{q) for the different spin-isospin modes, we will also find modes with huj{q) 
close to [{cqY + A^ p]^''^- These modes arise from the singularities in the form factors, 

i^.>,g)= ,5"7JT2T?\2 • (92) 

In the case of the pion, these singularities can be seen as the coupling of the pion 
to a heavier meson when the pion is spatially close to the nucleon, like in the elec- 
tromagnetic case where the pion couples to the photon via a p meson, see e.g. ref. 
1^. However, our model, with the monopole factors determined in the space-like 
sector {c\q\ > hu), is not appropriate to describe the physics near the singularities 
of the form factors where their true behavior may deviate substantially from the 
monopole form. To avoid this difficulty we exclude the "form factor" like modes 
{huo{q) ~ [{cqY + A^p]^/^) from the dispersion relations when calculating physical 
quantities, like cross sections and the A width, so that the form factors only contribute 
by their numerical values along the other spin-isospin modes. 

The short-range interaction contains the correlation parameters g'f^^, Qna^ 
g'/^j^, and the cut-off factor A^, and is an effective interaction that simulates more 
complicated interactions at short range, like exchange of heavier mesons and exchange 
of two or more mesons. Taking this interaction according to eq. ( |T3|) and excluding the 
form factor (A^ —>■ oo), this interaction becomes a 6 interaction in coordinate space, 
and is thus the simplest form of a short-range interaction. By including the form factor 
{Ag < oo), we take into account the finite size of the interacting particles. In reality, 
the effective short-range interaction may have some additional g-dependence via the 
(yf'-parameters, and may also depend on the nuclear density. However, in the present 
study we neglect such complications and take constant values for the (^'-parameters. 
The values of g'j^^^ dNA^ ^^11 known. Our preference is therefore 

to take Ag equal to either A^ or Ap, in order to reduce the number of free parameters, 
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and then vary the g' parameters to fit certain experimental data. For numerical 
reasons, it is more convenient to take = Ap than A^ = A^r- For A^ = A^ the 
singularity of Fg would lie close to the p-like branch in the spin-transverse channel, 
and because of the interaction between these two modes the p-meson mode would 
be substantially changed, and it would be difficult to exclude the "form factor" -like 
mode from the dispersion relations. To avoid these problems it is therefore convenient 
to take Kg = Ap. 



The value Ap = 1.5 GeV is a larger value than was used in ref. [|T^] were the A 
width in nuclear matter and its contribution from different decay modes was discussed. 
In that work A^ = 1.0 GeV was used. A g' set with a small cut-off factor should 
approximately correspond to a set with somewhat smaller g' values with a larger cut- 
off factor. This implies that the numerical values of the g' parameters in this work 



tot 



and ref. [T^] are only approximately similar. However, this does not affect how 



depends on the g' parameters, as was discussed in ref. [0. 

To determine the values of the g' parameters we calculate the p+p —>■ A+++n 
cross section in vacuum with our model. Keeping the previously discussed parame- 
ters fixed, we adjust g'j^^ to reproduce existing data. In fig. ^ we present calculations 
of da/dcos{6) at center-of-mass energy ^/s = 2.314 GeV together with experimen- 



tal data from ref. We find that we reproduce the angular dependence of the 

cross section well. In the calculations we have used the value g'^^ = 0.38. The cross 
section depends rather strongly on g'^/^, so therefore this value is quite well deter- 
mined in our model. In fig. ^ we present da/dt for the same parameter set, but 
with ^/s = 2.513 GeV, together with experimental data from ref. [^. We note that 



in ref. events with t or m equal to a given value were in the same bin. Thus, 



the experimental values have in fig. ^ been divided by a factor of two. Also at this 
energy the experimental data is reasonably well reproduced, except for the exper- 
imental peak value at t ~ —0.06 GeV^. In fig. ^ we present the calculated total 
cross section as a function of ^/s, together with two parameterizations and experi- 
mental data. The solid curve is our calculation, the dashed is the parameterization 



from VerWest-Arndt [23], and the dash-dotted line is a simple parameterization often 



used in BUU calculations p4|. The data points were estimated from fig. 2 of ref. 



23[] and can be found in references therein. As can be seen in fig. ^, we reproduce 
the energy dependence of the total cross section quite well, and substantially better 
than the simple parameterization from ref. At large ^/s we under predict the 

experimental points somewhat. Increasing g'^^ while keeping all other parameters 
fixed will increase the total cross section, but also flatten the angular distribution in 
da/dcos{6). Decreasing g'j^^ will reduce the total cross section and also flatten the 
angular distribution. 

The fit of the A^"*" cross section can be maintained if both f and g'j^^ simulta- 
neously are increased. For example, for = 7.2, /^^ = 12.2 and g'j^^ = 0.46 an 
equally good fit is obtained. However, the contributions from the spin-longitudinal 
and the spin-transverse channels are changed. For = 6.2 {g'j^^ = 0.38) these 
two channels contribute approximately equally, while at = 7.2 {g'jy^ = 0.46) the 
contribution to the cross section comes mainly from the transverse channel. 

The cross section cr(p + p A++ -|-n) does not depend on the correlation param- 
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eters g'p^j^ and (^aa; these remain undetermined by fitting cr(p + p — A++ + n). We 
can obtain some constraints on g'^^ and g'^^ by calculating the A width in nuclear 
matter at zero temperature, and identifying — r^*/2 with the imaginary part of the 
A-nucleus optical potential. In the A-hole model first developed for vr-nucleus scat- 
tering 1^ the authors used a A-nucleus optical potential with distinct contributions 



to the imaginary part, 

- 2 ImiVT) = - ^r^^"^^ - 2 Im(Kpread) • (93) 

The spreading potential was adjusted to fit the experimental results. For ^^C T^pread 
was found to be approximately 

^spread- [23 ±5-^(43 ±5)]^ McV, (94) 

Pn 

and rather independent of energy in the interval 100 MeV < T^^ < 250 MeV. Also 
in microscopic calculations of the A-nucleus optical potential |^ the authors have 



found that the imaginary part of the spreading potential is rather independent of 
energy, -40 < Im(V;pread) < -20 MeV, in the interval 50 MeV < pac < 400 MeV 
for a A on mass shell. Figure |^ shows the quantity [F^* — F^'^'^ + 5F^'^^'']/2, which 
corresponds to the spreading potential. The results for two g'^A values, 0.25 and 0.35 
are presented, using g'j^j^ = 0.9. We have performed the calculations presented in fig. 
^ at Pat = 0.75p^ to compare 0, ^ with the empirical values of ref. pSf. The A 



energy and momentum are estimated from the relations 



EA = ^^ + hu;M , pI = \p'f + q'. (95) 
5 irrij^ 5 



where q is the pion momentum determined from the pion kinetic energy T^- = hujT^{q) — 
mTrC^. The calculations shown in fig. ^ also include a binding correction of 20 MeV. 
The A width varies only slightly with (7^^, but depends quite strongly on g'^A^ ^^e 



ref. fl^. By comparison to the empirical points, it is seen that g^AA is quite well 
determined to be in the approximate interval 0.25 to 0.35. We have chosen to present 
our remaining results for the value g'AA = 0.35 

A value of (7^^ ~ 0.5-0.9 is often used in the literature. Some constraints on 
g'j^^j^ can be obtained from the low-energy (g ^ 0) Gamow- Teller response, as seen 
for example in (p,n) reactions. Several years ago there were reports that only 60% of 
the expected strength was found among the low-lying states. It was suggested that 
the low-energy strength is due to a strong couphng with AN~^ states, but the effect 
may equally well be explained by couplings to two-particle-two-hole states (see for 
example the review by Bertsch and Esbensen [^). A value of g'j^j^ ~ 0.9 is consistent 
with a rather weak coupling between low-lying Gamow- Teller modes and the AA^~^ 
states, while a smaller value, g'j^j^ ~ 0.6, leads to a renormalization by almost 40% of 
the low-lying GT strength due to the coupling to AN~^ states. The significance of 
the g' parameters is discussed further in ref. in connection with the calculation 
of the A width. 
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Since we perform our calculations at a constant density we will describe the nu- 
cleons by letting them propagate in a constant potential V^. In addition we take 
into account corrections in the nucleon energies, because of their interaction with the 
surrounding medium, by taking an effective nucleon mass as 



m 



N 



^1 

Po 



(96) 



This density dependent form is obtained in the extended Seyler-Blanchard model 
discussed in ref. p9|, ^ and leads to = 0.771 8m at at normal nuclear density, p = 
Po- That model employs an effective nucleon- nucleon interaction, with a Yukawa force 
modulated by a quadratic momentum dependence and an explicit density dependence, 
and is solved self-consistently within the Thomas-Fermi approximation. The model 
gives a good description of average properties of standard nuclear matter pO| . 

The A isobars propagate in the same manner in a constant A potential. For 
the modification of the real part of their energies, due to their interaction with the 
medium, we follow previous work on vr-nucleus scattering 
absorption , and nuclear response in the A region P3, P| 
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nuclear photo 
and take Va — Vn^ 25 
MeV at normal nuclear density. Note that in our formalism the potentials only enter 
as the difference V^ — Vn- This difference is found to be rather independent of the A 
energy and momentum p6[. The quantity V^ — Vn may be expected to vary somewhat 
with the nuclear density but in the present study we will use Va — V^v = 25 MeV also 
at twice normal nuclear density. 



4 Results and discussion 

In this section we present and discuss our results. For various specified densities 
and temperatures, we calculate the spin-isospin modes formed in symmetric nuclear 



matter. In sec. [4.1| we exhibit their dispersion relations, i. e. the energy-momentum 
relation uj{q), as well as their composition in terms of the unperturbed states. From 
these quantities, we calculate the width of the A isobar (sec. E]^) and the cross sections 



for collisions involving a A (in sec. |4.3| ). Furthermore, we discuss in sec. iA how the 



results could be incorporated in transport simulations of heavy-ion collisions, such as 



those carried out with the BUU model. In sec. we also discuss previous studies 
that have included some in-medium effects in transport calculations. 

The hadrons are confined within a periodic cubic box with side length L, and the 
coordinate system is aligned such that q, the momentum of the mode considered, is 
parallel to the z axis, q = qz. Since the box is finite, the momenta take on only 
discrete values, 

2nh . 27rh 

P = -J^J = -J^Ux, Jy,Jz) , (97) 

where j^, jy and are integers. The energy of an unperturbed nucleon- hole state is 
given by 

2 

^^^^-1 = Er,ip + qr) - E^(p) = -V + ^ , (98) 
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and for an unperturbed A- hole state (taking Fa = for simplicity), 

2 

Q QV 

hujAN-i = EA{p + q) - En{p) ~ \ + mA-mN ■ (99) 

Since the energies hu^N-^ and Tiojan-^ depend only on (and not and Py), it 
suffices to specify the quantum number jz in order to characterize a NN~^ or AA^~^ 
state. In reality the AA^~^ states have also a dependence, via the A width and the 
term p'^/2mA—p'^/2m*j^ which is neglected in eq. (^9]). We have taken this dependence 
into account approximately by making the replacement 

^ pl+<pl>ipz;T,pM) , (100) 

with the variance of the transverse momentum given by 

2 / m \ f d'^P± \ 2 I f d'^P± 

<P±> [Pz]T,pn) = J -^—^nN{p]T,pN) pj J -^—^nN{p]T,pN) , (101) 

where n^v is the nucleon occupation probability given in eq. (^Tf) . 

When calculating quantities like the energies and amplitudes of the spin-isospin 
modes, the total and partial A widths, and A cross sections, one can take the width 
Taat-i in the Lindhard function $a (see eq. ( |122| )) to be either identical to zero or 
equal to the total A width calculated self-consistently by an iterative procedure. In 
this section we will present results for both cases. We will refer to the former case as 
the reference case and to the latter as the self-consistent case. 



4.1 Spin-isospin modes 

From eq. ( |126| ) we calculate the energies of the spin-isospin modes that are formed in 
the interacting system, i.e. their dispersion relations. Figs. ^ and ^ display the real 
part of the dispersion relations for the self-consistent case at normal nuclear density, 
Pn = Pn = 0.153 fm~^, and zero temperature, T = MeV. In fig. ^ a number of 
different modes in the spin-longitudinal (vr-like) channel are apparent. Some of those 
are non- collective NN"^ modes (solid lines), which at zero temperature have their 
energies within the region 

2 

< Rehij < ^ + ^ , q<2pp , 
2m^ 

< Rehuj < -— ^ H , q>2pF . (102) 



2m*j^ m%j 2m% m%j 

Since we are presenting our results for a box normalization with a finite side length L, 
we obtain a discrete number of non-collective NN~^ modes. The total number of spin- 
isospin modes within the region ( |102| ) depends on L and tends towards a continuum in 
the limit L oo. In fig. ^ we show the NN~'^ component of the squared amplitude 
for all the modes at a fixed momentum, q = 300 MeV/c. It is clearly seen that all 
the squared amplitudes of the modes in the NN~^ region, < hcu < 185 MeV, are 
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dominated by a single NN~^ state (specified by Pz) and thus have a non-collective 
character. Another feature of the non-collective modes is that their energies are very 
close to the energies of the corresponding unperturbed states. 

Similarly, a number of non-collective AiV~^ states emerge in fig. which, for a 
fixed q and at zero temperature, have their energies constrained to a band, 

niA — rriN H < tie huj < niA — rriN H 1 • (lUo) 

"^A 2mA '^A 

The absence of collectivity is seen in fig. ^ where for all modes in the AA^~^ region, 
except a few in the upper half of the region, a single component (specified by pz) 
dominates each mode. The remaining modes show some collective behavior, but are 
mostly dominated by only two components. In fig. ^the solid curves show those AA^~^ 
modes that have a contribution larger than 0.5 from one or 0.9 from two components, 
while the remaining modes are represented by dot-dot-dot-dashed hues to indicate 
that they carry some collective strength. Also here the perturbed energies are close 
to the corresponding unperturbed ones. 

In addition, two collective modes are visible in fig. 0a. The lower one, represented 
by a dot-dashed line, starts at fiuj = m^c^ at g = and continues into the AN~^ 
region around q ~ 330 MeV/c. This mode is in the sequel be referred to as tti. The 
upper collective mode, displayed as the dot-dot-dashed curve, starts slightly above 
hu ~ mAC^ — m^vc^ at g = and approaches HcUt, = [(m^rC^)^ + (cg)^]'^/^ at large q. 
This mode is denoted n2. Fig. ||c gives an impression of the structure of the vfi and 
TT2 modes: the squared amplitudes of the pion and the sum of all NN~^ and AN~^ 
components are shown for a fixed momentum, q = 300 MeV/c. It is seen that the 
collective modes have contributions from all three types of interacting states, while the 
non-collective modes have only contribution from one type of state. Furthermore, we 
see from fig. |a and H& that the total NN~^ and AA^^^ contributions to the collective 
modes are made up from small contributions from all of the individual NN~'^ and 
AA^~^ states, respectively. 

Fig. ^ shows the squared amplitudes of the different components on the lower col- 
lective mode as a function of q. For small q the pion component dominates, while the 
AA^~^ component becomes dominant at around q ~ 330 MeV/c. At g ~ 500 MeV/c 
the lower collective mode has lost almost all of its collective character. 

In fig. 1^ we present the modes in the spin-transverse channel. The dispersion 
relations for the non-collective modes are very similar to the relations in fig. [7|a, while 
the collective p- meson like modes are different. There are two visible collective modes 
in fig. one dominated by the p-meson component is starting at hu = nipC^ at 
g = 0. This curve is represented by a dot-dashed curve and is denoted pi. The other 
one of A-hole character is denoted p2 and is displayed as a dot-dot-dashed line. This 
one starts slightly above hu k, rriAf? — mj^rc^ at g = and continues into the AA^~^ 
region around g ~ 200 MeV/c where it gradually looses its collective character. 

The dispersion relations presented in fig. |^ can also be calculated with Faat-i = 0, 
which gives purely real energies, Uu. The spin-isospin modes are then built up by 
stable A isobars having a fixed mass mA and so each mode will have a distinct 
energy. The self-consistent inclusion of the A width has several consequences. The 
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existence of a Breit-Wigner like distribution of A masses implies that the unperturbed 
AA^~^ energies no longer will be distinct, but also have a Breit-Wigner distribution 
with width T^^-i. This in turn leads to the same effect for the spin-isospin modes, 
where Re huu represents the centroid of the distribution and 2 Im hui, is its width. 
The A width also governs the decay rate of the isobar. Accordingly, the A isobars 
will have finite life times and so will the spin-isospin modes. The decay out of a mode 
will then occur by the process u —>■ AN~^, where the A will decay further. Thus Im 
hUi, contains information about the life time of the mode z/, before it disintegrates 
due to the decay of one of its constituent A components. 

The dispersion relations for the reference case, Faat-i = 0, are qualitatively similar 
to the real part of the eigenenergies obtained in the self-consistent case. Re hUi,, with 
differences mainly in the collective modes. When going from Faat-i = to the self- 
consistent treatment, the strength of the interaction between the states effectively 
weakens. This occurs because the inclusion of Faat-i = T^* causes the strength 
of the AN~^ states to be smeared out in the Breit-Wigner like distribution. As a 
consequence, for a fixed q, the energy of the pionic mode tti is somewhat raised and, 
by the same token, the energy of the upper mode 7^2 is lowered. The changes in energy 
is 0-50 MeV, depending on q and with the largest difference for large q. 

The imaginary parts of the eigenenergies hui, are presented in figs. ^ and 
The energies of the non-collective nucleon-hole modes are purely real, while the non- 
collective A-hole modes have imaginary parts that are close to half the corresponding 
widths Faat-i. For both collective modes, we find a large imaginary part, | Im huj^\, 
when the AN~^ components dominate the mode, and a smaller imaginary part when 
the mode is dominated by the meson component (compare fig. |^c and This 
implies that the collective modes will have longer life times when they are more 
meson-like than when they are AA^~^-like. 

In fig. ^ we present Re hui, for the self-consistent case at twice normal nuclear 
density at zero temperature, and in fig. ^ at normal density and temperature T = 
25 MeV.i Comparing the dispersion relations at normal (fig. ^ and double density 
(fig. P), we see that the main differences occur for the two pionic modes. The enhanced 
interaction at 2p% makes the two collective modes repel each other more strongly, 
which causes a lowering of the mode vfi and a concomitant rise in tt2- Another feature 
in fig. ^ is that the non-collective NN~'^ and AN~^ modes cover a larger region in 
the {u, q) plane since the nucleon chemical potential is larger at pjq = 2p^ than at 
Pn = P%- At T = 25 MeV (fig. |6) the collective modes tti and 7r2 are very similar 
to the modes at T = MeV (fig. |^, but there are many more non-collective modes. 
This is because the nucleon occupation probability, eq. (pTj), at finite temperatures 
allows the occupied nucleon-state to be found above the Fermi surface. 

^We limit our considerations of excited matter to rather moderate temperatures, because of the 
apphcations we have in mind. 
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4.2 A decay width 

The total A decay width in the nuclear medium, F^*, can be calculated according to 
eqs. (^Sp and (|TD|). In fig. |rU|we show the total A width for various nuclear densities 
and temperatures. The width is presented as a function of the invariant A mass, 
m = [E^ — (cp^)^]^/^/c^, for a fixed A momentum, pa = 300 MeV/c. Fig. [To| a 
displays Fa for the reference case, Faat-i = 0, while fig. |T0^ shows F^* when Faat-i 
is calculated self-consistently, Faat-i = F^*. 

The most noticeable difference between Fa* in the nuclear medium at normal 
nuclear density (solid curve) and in vacuum (dotted curve) is that in vacuum Fa* = 
Fa*^'' starts to increase from zero at the threshold m = mN + TriTr, while in the medium 
Fa* can be finite also for m < m^v + w^tt, since the A isobar can decay into a nucleon 
and a non-collective nucleon- hole mode, A —>■ N + NN~^. These modes have lower 
energies and the decay can thus occur below the free threshold. These effects are seen 



in fig. where also the partial contributions to the total A width are presented. We 
also see that at pn = p%, T = MeV, the total width in the medium. Fa*, is larger 
than Fa"^*^ up to m ~ 1400 MeV/c^. This effect is also mainly due to the new decay 
channels present in the medium, both A —>■ N + NN~'^ and also A ^ N + AN^'^ 
which begins to contribute at m ^ 1200 MeV/c^. 

Another feature seen, in comparing Fa'^'' with Fa* at pn = p%, is that Fa*^*^ 
increases more steeply than Fa*. This is an effect of the effective nucleon mass 
mlf < at Pat > 0. A A isobar with energy = [(mc^)^ -f- {cp^YY^'^ decays into 
a spin-isospin mode u and a nucleon N, and energy conservation yields E/^ = hu^+eN. 
By lowering m*^ the nucleon energy will be enhanced and thus the energy hui, 
of the mode u is reduced. This leads to a lower momentum of the mode u in the A 
decay, which in turn leads to a lower A width. 



The large enhancement of Fa* in fig. llOlat invariant masses up to m ^ 1250 MeV/c 
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for Pn = 2p^ occurs because the contributions from the nucleon- hole modes are pro- 
portional to the nuclear density, but it also refiects the fact that our parameter set 
leads to pion condensation at densities just above 2p%, as is manifested by the oc- 
currence of a spin-isospin mode lying in the region of nucleon-hole modes with a 



small real energy and a non- vanishing imaginary energy. As is also seen in fig. 11?^ 



almost all the contribution to Fa* comes from the nucleon-hole channels. The onset 
of pion condensation can be pushed up in density by increasing the values of the g' 
parameters g'j^^ and g'/^^, for example by making them density dependent. 

Increasing the temperature to T = 25 MeV has very little effect on the total A 
width. However, as seen in in fig. 0, the partial widths are considerably changed. In 
fig. |TT|a we see that at zero temperature there are mainly four contributions to the 
total width for m < 1400 MeV/c^: non-collective NN~^ and AA^^^ modes (long and 
short dashed curves, respectively) and the two collective modes vfi and 7^2 (dot-dashed 
and dot-dot-dashed curves, respectively). At low invariant mass the only energetically 
possible decay modes are the non-collective NN~^ modes. For m > 1100 MeV/c^ 
the A has enough energy to to decay to the mode tti and a nucleon above the Fermi 
surface. This becomes the dominating contribution to F^* from m ^ 1140 MeV/c^ 
up to m 1320 MeV/c^. The non-collective AA^~^ modes start to contribute at 
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m ^ 1200 MeV/c^ and they become dominant at m > 1320 MeV/c^. Accordingly, at 
zero temperature and normal nuclear density a A near resonance (m ^ m^) mainly 
decays into a nucleon and the pionic mode vfi, and the magnitude of the partial width 
is comparable to the free width. In addition, the partial width for decay into 
any of the non-collective nucleon-hole modes is approximately half of the free width, 
while the partial widths of the remaining decay channels are small at m ^ m^. At 
the temperature T = 25 MeV the situation is quite different. Concentrating still on 
a A near resonance, we find that the non-collective NN~'^ and AA^~^ modes give 
approximately equal and dominant contributions, while the contribution from the 
pionic mode is smaller. So, even though the probability for the A to decay is the 
same at T = MeV and T = 25 MeV, there will be very few decays to the pionic 
mode in the latter case. 

There are mainly two reasons for the reduction of at finite temperature. The 
first is that the mode tti loses its collective strength at a lower q value at finite 
temperatures, as compared to zero temperature. The second is that it is energetically 
possible for the A isobar to decay into non-collective AiV~^ modes at lower A energies 
at T = 25 MeV than at T = MeV. This is because there exits AA^~^ modes with 
lower energy at T = 25 MeV than at T = MeV, and because the nucleon formed 
in the A decay can have energy less than Ci;' at finite temperatures. With several 
more decay channels contributing to F^* at T = 25 MeV the contribution from 
is reduced, since the total width is almost unaffected by the change in temperature. 

The decomposition into partial widths corresponding to specific spin-isospin modes 
z/ is unambiguous when F^^-i vanishes. However, in some cases it is not possible to 
determine uniquely the nature of a specific mode u. One such case is when the collec- 
tive mode u enters the non-collective AA^~^ region and loses its collective strength. 
This transition is of course gradual and there is some arbitrariness involved in de- 
termining when the mode is no longer collective. This gives rise to a corresponding 
(small) uncertainty in the value of the partial widths and . Furthermore, 

the interaction between the pionic mode and the non-collective AA^~^ modes gives 
rise to some collective strength among a few of the AA^~^ modes for certain small 
intervals of the momentum q, see fig. ^. It may be debated whether the contribution 
from these collective modes to F^* should be associated with F^^ , V^, or some- 
thing else. In figs. ^ and we have indicated the contribution from all such cases 
by error bars on the partial widths. 

From the spin transverse channel we get a contribution to F^* from the non- 
collective NN~^ and AA^^^ modes, and two collective modes pi and p2- However in 
the range of invariant masses presented in figs. |ll] and |12] the contributions F^^ and 
F^^ are negligible or small and therefore not displayed in the figures. The partial 
quantities F^^ and F^^ contain contributions from both spin longitudinal and 
transverse channels. 



In fig. we present the partial widths based on eq. ([73| ) with the approximation 
([riD, where the A width is included self-consistently. In fig. |l2|a the partial widths 
are shown at zero temperature and normal nuclear density. Comparing the widths 
for a decaying A near resonance (m fa m^) with the reference case F^^^tv-i = in 
fig. |TT|a, we see that F^^ is substantially reduced, F^^ is rather unaffected, while 
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we get a substantial contribution from the partial widths. The relatively large 
value of already at m ^ m/^ occurs because it is now energetically possible to 
form the mode 7f2; since it may exist at lower energies, as compared to the case with 
Faat-i = 0. The reduction of F^^ is understood from the facts that the total width 
is rather unaffected by the inclusion of = F^* and that new competing decay 

channels have opened up. 



4.3 A cross sections 

In this section we present results for cross sections for the processes 

N + N ^ ^ + N , N + ^-^ N + N (104) 

and 

j = l,2, (105) 

where iXj denotes the collective spin-isospin mode in the spin-longitudinal (pion like) 
channel. 

In vacuum the cross sections for the two first processes will depend on the trans- 
ferred energy and momentum, {uj,q), in the relativistically invariant form uj'^ — q^. 
It is therefore convenient to perform the calculations in the center-of-mass system of 
the colliding particles, where the total cross section only will depend on the total en- 
ergy, y/s, of the colliding particles. In the medium, however, the spin-isospin modes, 
and thus the effective interaction, no longer depend on (cj, q) in the invariant form 



a;^ — q^, see fig. 0. This implies that the total cross sections, apart from depending on 
the total center-of-mass energy, will have additional dependences on the momenta of 
the colliding particles. However, since our formalism for calculating the isospin modes 
and the effective interaction is non-relativistic, the calculations should be performed 
in the rest frame of the medium (where the spin-isospin modes are calculated). In or- 
der to study the medium effects, we wish to compare to experimental and calculated 
cross sections in vacuum. In this section we therefore calculate and present all cross 
sections for the special case that the center-of-mass system of the colliding particles 
is identical to the rest frame of the medium. 

In fig. Owe present the differential cross section da/dcos{6cm) for the process 



iV ^ A + A^ at a center-of-mass energy of ^/s = 2314 MeV. The figure shows different 
combinations of density and temperature. The calculations of fig. correspond to 
Taat-i = 0, while the results of fig. |13|& represent the self-consistent case. Comparing 
the results in the nuclear medium with the results in vacuum (dotted curve), we 
find that in normal nuclear matter at zero temperature, fig. 0, the cross section 
is somewhat enhanced at forward and backward angles and slightly suppressed at 
cos(^cm) ~ 0. The enhancement at forward and backward angles occurs because the 
pionic mode tti is lowered in the medium, see fig. |^. The scattering matrix contains 
terms proportional to [hu — hu!^{q)]~^ for each spin-isospin mode, see eq. (Q). The 
dominant term in the sum over the modes u is the pionic mode u = fci, and its 
energy huji{q), keeping q fixed, decreases when the nuclear density increases. Hence 
the difference huj — huji{q) becomes smaller at forward or backward angles, and the 
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cross section is enhanced in the medium. This effect is also responsible for the large 
enhancement of the cross section at twice normal nuclear density seen in fig. At 
the temperature T = 25 MeV the differential cross section is only slightly reduced as 
compared to zero temperature. 

Including the A width self-consistently, via F/^tv-i in eq. ( |122[ ), reduces the cross 
section only slightly at normal density, fig. 13b, but to a significant degree at twice 
normal density, fig. |TB|6. This is because the pionic mode is regularized by the imag- 
inary part of the pion optical potential. As a consequence, the pionic mode has no 
longer a single well-defined energy, but instead there is a Breit-Wigner-like distri- 
bution of possible energies centered around Re huji with a width 2Im hui. In the 
former case, when the pionic mode had a well-defined real energy, ui, the transferred 
energy and momentum comes close to this energy, which causes a large enhancement 
of the cross section, but in the latter case the transferred energy and momentum will 
only be close to some of the possible energies around Re huji and thus only pick up 
a fraction of the total strength in the pionic mode. Technically, this is taken care 
of by the imaginary part of hui which increases the difference huj — huji{q) in the 
self-consistent case. The reason why the effect of including FAiv-i = T^* is so much 
larger at pjv = 2p^ than at pat = is mainly that the transferred energy oj is closer 
to Re uji at 2p^, but also that Im uji is larger at 2p^. The magnitude of Im ui 
depends on the magnitude of F^* and at 2p^ the width F^* is larger than at p^. 

In fig. [1^ we present the total cross section for the process N + N — A + 
as a function of the center-of-mass energy ^/s of the two colliding nucleons. The 
different curves correspond to different densities and temperatures. Fig. fL^ i displays 
(Ttot for Faat-i = 0, while fig. |l4b shows atot for the case when Faat-i is included self- 



consistently. In fig. |T^ we notice that the total cross section is enhanced by a factor 
2-3 at twice normal nuclear density, as compared to the value in vacuum, while at 
normal nuclear density the cross section is similar to the vacuum value with only minor 
changes. As for the differential cross section, the enhancement of atot at pat = 2p^ 
originates in the softening of the pionic mode. The self-consistent inclusion of Faat-i 
has two major effects. At normal nuclear density the cross section is enhanced in 
the threshold region {y/s ^2.1 GeV) and somewhat reduced at energies ^/s ^2.2- 
2.4 GeV. This effect arises from the possibility of creating A isobars with masses lower 
than mAT -l- m^^ in the nuclear medium and it is technically included by taking the A 
width in eq. ( P5| ) from the nuclear matter calculations presented in sec. [4.2| . At twice 
normal density the cross section is reduced at energies a/s i=a 2.2-2.4 GeV by about 
10 mb, as compared to the when case F^tv-i = 0, while it is enhanced significantly in 
the threshold region. The origin of the reduction has already been discussed above 
in connection with the discussion of the differential cross section. The enhancement 
at low a/s occurs because the A width is very large at low A energies (see fig. p!0|). 
The increase of the A width at low A masses and = 2p^ is partly an effect of the 
system being close to the onset of pion condensation. The associated enhancement of 
(Ttot is consistent with the picture that when the system is close to pion condensation 
the A isobars having low energy quickly decay into low-energy pionic modes in the 
NN~^ region. The effects at 2p^ may be quantitatively changed somewhat by taking 
into account Re Sa and its density dependence. 
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The total cross section for the reverse process A^ + A N + N is presented in fig. 
TBI for a A with mass m = 1230 MeV/c^. As in fig. |14|, part (a) presents the results 



for r^^-i = and part (b) corresponds to Faat-i = T^*. Similar features as in fig. [T| 



are noted in fig. |T5|. The vacuum cross sections are only slightly reduced at pn = p%, 
but for F/N^Tv-i = greatly enhanced at pn = 2p^ while F^^-i = F^* leads to a 
cross section of similar magnitude at all presented densities for m = 1230 MeV. 
Furthermore, we note that (7{N + A —>■ N + N) is singular at threshold. This 
singularity originates from the factor \pi\/P in the cross section, eq. (0), which in 
the A^+A center-of-mass system becomes proportional to l/|Pi| when i/i — > mi+m2- 
As this limit is approached, \pi\ will tend to zero and hence the cross section grows 
infinite. Finally, we note that the cross sections at A masses m < 1230 MeV/c^ are 
smaller than at m = 1230 MeV/c^. 

In fig. |l^ we present, for various nuclear densities and temperatures, the cross 
section for the process p+7rt — > A++, when Faat-i is included self consistently. In 



figs. p!6| a and p!q & we present the results for the lower pionic mode, vfi, and in figs. [T( 
and |TB|(i for the upper pionic mode, The factor p^^m'^), defined in eq. (^), is in 
figs. |l^a and |l6]c calculated from the free A width, while the total width in nuclear 
matter has been used in figs. ^6 and |TB|d. 

The cross sections are presented as a function of the invariant A mass, m, which 
is equal to the total energy in the p+tTj center-of-mass system. The cross sections 
show a strong resonance peak at m for all densities and temperatures. This 

shape is in eq. (^) determined by the factor p3(m^). The magnitude of the cross 
section is determined by this factor and the factors |/i''*(AA^, a;j,)p and frci/c. In 
vacuum, where only the real pion dispersion relation exists in the spin longitudinal 
channel, the factor is proportional to {cqY/hu.„{q). At finite densities the pion 
contribution to will be mixed with contributions from AA^^^ and NN~'^ states 
which depend very weakly on the momentum and are approximately proportional 
to g'^/^PN and g'j^^pN, respectively. On the lower pionic mode the pion component 
dominates for low q values, while the AA^~^ components dominate at larger q, see fig. 

Therefore, the cross section for the lower pionic mode at finite densities will be 
similar to the vacuum cross section for small values of q or m, while it will be smaller 
than the vacuum value for larger q. This is seen in fig. p!6|a, where the cross sections 
at finite densities deviates from the vacuum values for m > 1200 MeV/c^. At normal 
nuclear density and zero temperature the momentum q is approximately 270 MeV/c 
at the resonance peak, m ^ m^. At this momentum the AA^^^ component in the 
pionic mode is quite substantial, see fig. At twice normal density the contribution 
from the AA^~^ component is approximately doubled, which is the main reason for 
the decrease of the cross section at the resonance peak. Moreover, the value of the 
relative velocity will be changed at finite densities. There are two competing effects 
on the lower pionic mode. The velocity of the pionic mode is lower in the medium, as 
compared to a pion in vacuum, since hui increases less steeply with q than hu.,,. On 
the other hand, the velocity of the nucleon will be larger, since m = m/\ = eN{q) + huJi 
is obtained at a larger value of g, because hCoi is lowered in the medium. The net 
results are seen in figure |1^. At finite temperatures the cross section is similar to the 
case of zero temperature. 
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In fig. we have used the A width in the nuclear medium to calculate the factor 
P3(m^), rather than the free width, used in fig. |l6|a. In the nuclear medium the A 
width is larger than the free width up to invariant masses around 1400 MeV/c^. This 
explains the reduction of all cross sections at the resonance peak in fig. ^6 relative 
to fig. p!6|a, since at the resonance peak the factor P3{m\) is proportional to I/Fa- 
Note also that the cross sections for the lower pionic mode in fig. [T6| 6 are substantially 
enhanced for low invariant masses. Also here the effects at 2p^ may be quantitatively 
changed somewhat by taking into account Re Ea and its density dependence. 

On the upper pionic mode the relation of the strength of the pionic and AA^~^ 
components are opposite, and thus the cross section will approach the vacuum value 
at large q, while it will be smaller at low q, see figs. |T^ c and [T6|(i. Note also that the 
invariant A mass will always be larger than on this mode, since huj2 starts just 
above m^/c^ — m^v/c^. In the limit g — the cross section diverges, since in this 
limit also the relative velocity approaches zero. 

4.4 Implications for transport descriptions 

In this section we discuss how our results could be incorporated in a dynamical trans- 
port simulation of a heavy-ion collision, and some of the consequences this could lead 
to. We will carry out this discussion within the framework of a standard quasipar- 
ticle description that propagates nucleons (A^), delta isobars (A), and pions (tt) as 
ingredients and we will refer to this as the "standard" transport description. In such 
a model, many of the A and tt properties, such as decay widths, cross sections, and 
dispersion relations are usually taken as the properties in vacuum. In this section we 
will discuss how the vacuum properties can be replaced by in-medium properties in 
a consistent way. 

In a transport simulation of heavy-ion collisions, a nucleon-hole {NN~^) excitation 
is produced by promoting a nucleon from below to above the Fermi surface. This can 
occur as a result of the nucleon colliding with another particle. The new NN-'^ state 
is then by construction non- collective and unperturbed, i.e. its energy is given by 
-^particle — -E'hoic- Thus the non-collcctive NN^^ spin-isospin modes that we have found 
in nuclear matter are already incorporated in the standard transport description. But 
the energy of the NN-^ state is given by the respective quasiparticle energies and so 
it is not quite correct and should in principle be slightly changed, in accordance with 



the energies of the non- collective NN~ modes presented in sec. ^Aj. However, the 
energies of the non-collective NN~^ states are only slightly shifted from the respective 
unperturbed energies, and it should thus be a good approximation to neglect this 



change of the energy. Similarly, the non-collective AA^ ^ modes in sec. ff^ correspond 



in a transport description to the conversion of an individual nucleon to a A isobar. 

The incorporation of the two collective spin-isospin modes is more involved. These 
modes can be regarded as separate particles of pionic character, vfiand 7i2, and treated 
in a manner analogous to the standard treatment of the pion. Since the pion is then 
fully included in the description, it should no longer be treated explicitly. The prop- 
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agation of the two collective pionic modes is governed by the effective Hamiltonians 



Hi{r,q) = Rehui{q;pN{r),T{r)) = hui , 

H2{r,q) = Rehu2iq;pNir),T{r)) = ^2 , (106) 

where fiuji and huj2 are the energy- momentum relations for the lower and upper 
collective modes discussed in sec. ^]l] and displayed in fig. ^ for p^r = and T = 0. 



Note that the spatial dependence of Hi{r,q) is incorporated by representing pwir) 
and T(r) as local quantities. Moreover, in the collision term the process for the 
production and absorption of pions in the standard description, A ^ + tt, should 
be replaced by the two distinct processes 

A ^ iV + 7fi and A ^ + tts • (107) 

The A decay is governed by the A decay width in the medium to these two specific 



channels, . These partial widths, presented in sec. [4.2| , should be employed in the 
same manner as the free width, i.e. they describe the probability for the A isobar to 
decay into a nucleon and a pion. The only difference is that several collective pionic 



modes are available in the final state. The reverse processes in ( |107| ) are characterized 



by the cross sections that were discussed in sec. ^73 . 

The self-consistent inclusion of the A width in the calculation of the spin-isospin 
modes encompasses decay processes like 

However, since such processes are already explicitly contained in the transport simu- 
lation by processes like 

it would not be correct to include the entire self-consistent A width when calculating 
the collective modes to be used in the transport description. Instead it is more correct 
to use the results obtained with Faat-i = 0, both for the energies of the modes, hojj, 
and for the partial A widths to be used in the decays A N+ttj. Still, it is important 
to use the self-consistent width when calculating the cross sections for processes like 
N + N — > A + A^, where the spin-isospin modes provide an intermediate effective 
interaction. 

Although the collective pionic modes can thus be effectively treated as ordinary 
particles, the fact that their wave functions contain components from vr, NN'-^ and 
AA^~^ states makes it difficult to picture them in a physically simple manner. Fortu- 
nately, their specific structure is irrelevant, as as long as these quasiparticles remain 
well inside the nuclear medium. First when such a quasiparticle penetrates a nuclear 
surface and emerges as a free particle is it physically meaningful to determine what 
kind of real particle it is. The gradual transformation of the collective quasiparticle is 
automatically taken care of within the formalism, because as the density is lowered, 
Pn 0, the pionic modes will acquire 100% of either the pion component or the 
AA^~^ component, depending on u and q. That is to say, they will turn into either 
a free pion or an unperturbed AA^~^ state. There remains the practical problem of 
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how to represent an unperturbed AA^~^ state when 0. However, we anticipate 

that only a very small fraction of the pionic modes will emerge as unperturbed AA^~^ 
states. For the lower pionic mode, it is only for large q that the AA^~^ component 
will dominate at low densities, but at large q the lower pionic mode starts gradually 
to lose its collective character and there is therefore a low probability for creating the 



lower pionic mode at large values of g, see fig. |T^. On the other hand, for the upper 
mode it is at low q that the A ^-component will dominate. In cold nuclear matter 
the A decay at low energies (which corresponds to low q) is strongly reduced, because 
the nucleon produced in the decay is Pauli blocked. In an actual heavy-ion simulation 
this reduction is smaller, but we still expect that the number of upper pionic modes 
at low q values will be quite small. These effects will be further investigated and 
reported in a subsequent paper. 

As was discussed in sec. the partial decay width is reduced at T = 25 MeV 
compared to T = 0. Thus in a heavy- ion collision there will be very few pionic modes 
TTi created in the hot region, while in colder regions the A isobars near resonance 
decay mainly into the pionic modes. In sec. |4.2| we also found, at twice normal 
nuclear density, a large enhancement in the A width at m ^ 1200 MeV, signaling the 
onset of pion condensation. Although our parameter set leads to pion condensation 
at > 2p^ for an equilibrated infinite system, it is not clear that the effect will be 
seen in a nuclear collision, since the region where such high densities may be created 
is rather small, and exists for a only short time. These effects will also be further 
investigated. 

The calculations of Fa presented in figs. pUf^T^ take account of the Pauli blocking 
of the nucleon in the A decay A ^ N + v. In a transport description the Pauli 
blocking of the nucleon is treated explicitly and should thus not be included in the 



width of the A. In fig. 17 we present the total and partial widths, without Pauli 



blocking of the nucleon, at different temperatures for the reference case (Faat-i = 0). 

The total A decay width, has apart from the partial contributions F^% also the 
partial contributions F^^ and F^^ . The partial width F^^ gives the probabil- 
ity for the A to decay into a nucleon and a NN~^ state. In a transport description, 
this implies that we initially have a A and after the decay process we have two nu- 
cleons above the Fermi surface and a hole left in the Fermi sea. But this is the same 
process as if the A would collide with a nucleon below the Fermi surface to give two 
nucleons above the Fermi surface. This process is normally already included in the 
collision term in a standard transport description, and the probability for such a col- 
lision is given by the cross section for the process IS. + N ^ N + N . In a transport 
description it is therefore not correct to both include a A decay according to F^^ 
and a collision term with A + — > + A^. Instead, the correct procedure should 
be to exclude F^^ and modify the cross section (t(A -|- A^ ^ A^ + A^) to be the 
in-medium cross section. Calculations of such in-medium cross sections is discussed 
in sec. In the same way, F^^ should be excluded in a transport description, 
and cr(A + A^ ^ A + A^) be the in-medium cross section. 

As mentioned in the introduction, the in-medium effects on the cross section for 
N + N A + N have previously been investigated in ref . [O] . That work corresponds 
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to zero temperature, T = 0, and taking T^jy-i = in $a, eq. ( |122| ). Our results in 



figs. and \L% i are in qualitative agreement with those of ref. |T^, i.e. a large 



enhancement of the cross section at high nuclear densities caused by a lowering of 



the pionic mode. However, as was also pointed out but not pursued in ref. |12], the 



imaginary part of the pion optical potential regularizes this effect. This is seen by 



comparing parts a and b of figs. O and 14, where it is shown that the cross section at 



twice normal nuclear density is considerably different if the A width is included self- 
consistently in the formalism. This points towards the importance of incorporating 
the A width consistently within the model for obtaining the correct magnitude of A 
cross sections at large nuclear densities. 

Since the N + A N + N cross section is not experimentally known one has 
traditionally in transport descriptions used the N + N A + vacuum cross section 
and detailed balance. As pointed out by previous authors, for example ^ it is 



then important to take into account the finite width of the A, i.e. use an expression 
according to 

da{N + N + N) I 1 da{N + N ^ A + N) 



dVL Nipf p'iim?) dVtdm? 
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where p\ and pf are magnitudes of the initial and final momenta in the center-of-mass 
system of the A^A system, A^f is a spin-isospin factor and p^{m?) is defined in eq. 
(pSl). The condition in eq. ( |108| ) is automatically satisfied within our formalism. 

In refs. []T3|, Q the authors included some medium effects of the pionic modes. In 
these works a simple form of the pion polarization function were used. This simple 
form originates from approximating the continuum of non-interacting AA^~^ states to 
a single state with energy hu/^ = m/^ — m^r + (gc)^/2mA, and neglecting the nucleon- 
hole states. When the interaction between the AA^~^ state and the pion is turned 
on two collective states emerge (but contrary to our model, no non-collective AA^~^ 
states are left). This approximation leads to some differences as compared to the 
more complete treatment in this paper. As seen in fig. 0, the lower collective mode vfi 
disappears {i.e. loses its collective strength) when it enters the non-collective AA^^^ 
region. This is not the case for the approximation used in [|13l n, where the lower 



pionic mode exists for all momenta q. Furthermore, the change in the dispersion 
relations hujp{q), relative to the relations in vacuum, are somewhat overestimated 
when the approximate form of the polarization function Tl^^ is used, as compared the 
the more complete treatment in this work. Finally, no justification was given for the 
omission of the A width in the pion polarization function. 

To avoid the problem of how to treat a non-interacting AA^~^ state penetrating 
the nuclear surface {i.e. when 0) in the transport description, the authors in 



1311 derive effective dispersion relations corresponding to asymptotically free pions 
or AA^~^ states. However, the energies in the effective dispersion relations are quite 
different from the original ones, especially for large momenta q. At g ~ 700 MeV/c 
the difference is as large as 150-200 MeV. 



In ref. [|T3[ the authors propagate only the collective mode that corresponds to 
the asymptotically free pions (71^2). The other collective mode is identified with the 
propagation of individual A isobars. This appears to be incorrect, since (as we have 
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argued) the propagation of individual isobars should be identified with the (remain- 
ing) non- collective AA^~^ states, and both collective modes should be treated on an 
equal footing in the transport description. The authors in ref. also used the en- 
ergy of the collective mode that corresponds to the asymptotically free AA^~^ states 
{712) to estimate the density dependence of the A potential. This had a large effect 
on their results. However, our calculations show that the non-collective AA^^^ modes 
have their energies very close to the unperturbed energies, and thus change very little 
with the nuclear density. We therefore feel that it is inappropriate to use the proper- 
ties of the collective modes to deduce any density dependence of the A potential to 
be used for the explicit propagation of uncoupled A isobars. 

The authors of ref. [14| do treat the two collective modes on an equal footing, prop- 
agating both of them as quasipions in accordance with the dispersion relation implied 
by their model, and their model also contains explicit propagation of uncoupled A 
isobars. They found that only a small fraction of the quasipions approach free A-hole 
states as they penetrate the nuclear surface. These quasipions were approximated by 
on-shell pions at the surface, by changing the momenta of the quasipions. 

Furthermore, in ref. ||14[ the authors calculated partial widths for the A decay to 
a nucleon and a quasipion. The authors state that there are also contributions to the 
total A width from decay processes as A — iV -|- NN'-^ and A — -|- AA^~^, which 
are not taken into account since these processes are already included in the transport 
simulations. However, by using the simple form of the pion polarization function, 
the AA^"^ continuum is compressed to a single state. But this state contains all the 
strength of the AA^~^ continuum, and therefore the partial A width for this channel 
will also contain some contribution from the A — -|- AA^~^ decay. 

Neither of the works |T3[ nor [|I^ have taken into account any modifications of 
cross sections for processes involving a A. 



5 Summary 

We have investigated the properties of spin-isospin modes in an infinite system of 
interacting nucleons, A isobars, and vr and p mesons at various densities and temper- 
atures. The aim has been to derive and discuss quantities that can be incorporated 
in a transport description of a heavy-ion collision, by use of a local density and tem- 
perature approximation. 

Within the random-phase approximation we have derived dispersion relations for 
the spin-isospin modes and the amplitudes of the their components. While the dis- 
persion relations yield the energy-momentum relation of each mode, the character of 
the modes is determined by the amplitudes of the different components. In both the 
spin-longitudinal and spin-transverse channels, we find two collective modes, while 
the remaining modes are non-collective in their nature (except in limited regions in q 
for a few specific modes). The non-collective modes correspond in a transport descrip- 
tion to propagation of uncoupled nucleons and A isobars, while the collective modes 
correspond to propagation of quasimesons. These quasimesons can be incorporated 
in a transport description in a manner analogous to how real pions have been incor- 
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porated in standard treatments based on vacuum properties. One notable feature of 
the lower pionic mode in the spin- longitudinal channel is that it gradually loses its 
collective character when it enters the region of of non-collective modes. This implies 
that a A isobar with sufficiently high energy cannot decay to this mode. Previous 
works [0, |T^ that included some in-medium effects employed a simpler model for 
the collective modes. In that simpler model the lower collective mode exists for all 
momenta q and can thus be excited by a decaying A at any energy, in contrast to our 
more refined results. 

The decay of a A isobar into a nucleon and a spin-isospin mode is governed by 
the partial A decay widths. Therefore, we have calculated total and partial A decay 
widths within the model. At twice normal density the total A width is significantly 
enhanced at low A energies. However, this enhancement is mainly associated with 
the decay to non- collective nucleon-hole modes, while instead the partial width for 
the decay to the lower pionic mode vfi is reduced. This effect is also different from 



the previous works |T3|, |14[ where the nucleon-hole channel was neglected in the 
calculation of Fa. At finite temperatures up to T = 25 MeV the total A width is 
almost unaffected, while the partial widths change somewhat with T. The dependence 
is stronger when Faat-i = than for the self-consistent case, Faat-i = F^*. 

The partial A widths representing decay to non-collective modes correspond in 
transport models to processes like A + N + N. Since these processes are already 
explicitly included in the transport description, these partial widths should be ignored, 
while the corresponding cross sections should contain the in-medium modifications. 
Examples of such in-medium cross sections have been presented. At center-of-mass 
energies y/s ~ 2.3 GeV, we have for the reference case F^at-i = found a large 
increase in cr(A^ + A^ A + N) at = 2p% as compared with p^r = p%, in agreement 
with ref. |T^. However, this effect is substantially reduced in the self-consistent 
treatment, Faa^-i = F^*. Instead, the cross section is significantly enhanced at low 
center-of-mass energies, {^/s ~ 2.0 — 2.2 GeV). This was not investigated in ref. 



T^ . Furthermore, the cross section of the process A-|-A^— >-A^ + A^is reduced when 
rA7v-i is included self-consistently, compared to when Faat-i is zero. We have found 
that all the calculated cross sections are almost independent of temperature up to 
T = 25 MeV. 

In a forthcoming paper we will incorporate the in-medium effects presented in this 
paper into a microscopic transport model. We will study the importance of the these 
in-medium properties in heavy-ion collisions and compare to the previous treatments 



T^, [TB|, We expect that the A production will be enhanced, especially at low 
center-of-mass energies of the colliding nucleons, but we also expect that the A decay 
to a pionic mode will be somewhat reduced. This could lead to a faster thermalization 
of the system and possibly slightly less real pions produced. The net effect, however, 
is difficult to predict without an explicit transport simulation, since an average over 
different densities and temperatures will be taken, and the pionic modes that will 
escape the system as real pions will be produced mainly at the nuclear surface where 
the in-medium effects are small. 

The model presented in this paper constitutes a more consistent way of obtaining 
and incorporating in-medium effects in transport descriptions of heavy-ion collisions 
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than previous works [0, 13, Ti]. However, also in the model presented in this paper 
some approximations and assumptions have been made, and there is room for fur- 
ther improvements. We expect that our model will be applicable in transport models 
up to moderately large bombarding energies of about 1 GeV per nucleon, since it is 
assumed within that the density of A isobars and pions is relatively small. Further- 
more, we have so far not included any density dependence of the coupling constants, 
g' correlation parameters, form factors, and meson masses, such as may result from 
a possible partial chiral restoration, since we feel that such effects are not very well 
known, at the present stage. 

This work was supported by the Swedish Natural Science Research Council, by 
the National Institute for Nuclear Theory at the University of Washington in Seattle, 
and by the Director, Office of Energy Research, Office of High Energy and Nuclear 
Physics, Nuclear Physics Division of the U.S. Department of Energy under Contract 
No. DE-AC03-76SF00098. 

A Solution to the RPA equations for spin-isospin 
interaction 

In this section we show how to obtain the eigenenergies of the RPA equations (^) in 
the spin longitudinal channel using the interactions defined by eqs. (^) to (|1^). The 
eigenenergies in the spin transverse channel is obtained analogously. 

We wish to find an spin-isospin excitation propagating with momentum q and 
isospin A. To such a spin-isospin mode that has a momentum q, only the baryon 
pairs that has the relative momentum q will contribute. For this purpose we need 



to restrict the summation over all baryon states in eqs. (^) and (|36D to those that 
will have the relative momentum q. In the same way we restrict the sum over meson 
states to those with momentum q and isospin A. Therefore we take 

Xjfc Xjk{uJ,q,X) 6p^^p^+q , (109) 
Zr -> Z{uj,q,\) Sq^^q Sx,,x , (110) 

Wr ^ iy(^,g,A) 5g^„_g 5a„,-A . (Ill) 

In q^j^ (eq. (H)) we will neglect the term pj^ {Tiuj) / {rriNC^ + Tioj) which is small in 
NN~^ or AA^~^ states, since the hole momentum is small, and take 

^2 



rriNC 



(112) 



Furthermore we take 

^/s = (^Ar(Piv) + kuj)'^ - {cq + cpj^f ^ [niNC^ + Re hujf - c{qf = ^il . (113) 
With these approximations we can make the ansatz, 
V ( \\ f^^Y^^ x{tj,tk; q; \, uj) 

Xjkiu, q,X) = l — \ — ■ — . ^ , ^ T— i^jkiq, -A , 114 

\LJ et{Pk + q)-et^{pj^) + Ejk-hu; 
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We consider the case when a;(3/2, 3/2; qr; A, lo) = 0. We then write x(l/2, 1/2) = , 
x(l/2,3/2) = a;(3/2, 1/2) = x^, and analogously for other quantities. The RPA 
equations can then be written in the matrix form 



1 + }^^;NNJ^N^N yyiVA^y^A^A 



TV 



and 



Z 
W 



hUJTT + ^ 



X 

x' 







<A„,A^A_A 



with 



fNaf. 



NI3 



, ,\ • DO 77' 



2m 



where a,j3 — N,A, and where we have defined the Lindhard functions 

n{p) - n{p + q) 



(115) 

(116) 
(117) 

(118) 
(119) 

(120) 



p + 9)^/2m;^ - p'^/2m*j^ - hw 
$^(u;,qf) = $(-,-;u;,qr) + $(-,-;u;,qr)=(^-^j X] 



(121) 
(122) 



with 



^^AN — 



Se 



- ^ + Am + l^/(£^,pi) - ^T^^-^{et,pi)/2 - fiw (123) 

{Pa? 



P^ 



+ Am + Vf{e^,pl) - iT^N-i{s^,pl)/2 + huj (124) 



AN 



2mA 2m 



N 



P^ 



Am = mA - mAr; £s = "^at + ±hu;; p^^p± q . (125) 

The numerical factors M.^ — A and A^'^ = 16/9 originates from the spin-isospin 
summation. 

The eigenenergies are obtained from 



yyANj^N^N 



det 



1 + W^^M^^^ 



0, 



(126) 
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with the approximation T^^-^i^h + ^t^) ~ ^AN-^i^h + Re hu). The eigenenergies of 
the auxihary equations are obtained from 



det 



1 + y\;^Nj^N^N -^iVA^y^A^A 
Y^AN^N^N I ^ yyAA_y^A^A 







(127) 



Since the normal RPA equations (|38| ) and the auxiliary RPA equations differ only in 



the matrix A^^\ we can use the relationships $ 



N, 



and $A(c(j*) = (^^{uj)* to show that if is a solution of eq. 
solution of eq. ( |127| ). 

From the normalization condition (^) we obtain 

±1 = M^7]nx*^xn + M^7]Ax\xA + Z*Z -W*W 



, then = cu* is a 



x*j^Xn \M r]N + 
1 



1 + y^NN^N^^N 

yyNAj\yiA^A 

1 



1 + W^^A<^<I>^ ^' 



:i28) 



with 



VA{uJ,q) 



he 



n{p) - n{p + q) 



d 



P i 2m* 



2m*. 



dhuj 



^N{uj,q) 



(129) 



he 



E 
p 



n{p)[l + idTAN-ii£h+ Re huj^p^) /2dhLo\ 



n{p) [1 + iSFAAT-i {sh — Re huj, p^)/2dhuS\ 



[6e 



AN\ 



d 



dhu 



$A(^,g) 



(130) 



where dV^ /dhuj ^ 0. 

The quantities of the auxiliary equations are related to the ordinary quantities by 



Xh(UJ 



Xh(UJ] 



b = N,A; 



Ziuj* 



-Ziuj) 



-W{uj) . 

(131) 

It is straightforward to show that the solutions defined by eqs. (|114|) , ( [116| ), ( |117| ), 
( |126| ) and (|128|) satisfies eqs. (|38|) , and that the corresponding auxiliary solutions 
solves the auxiliary equations. 



B The effective spin-isospin interaction 

A spin-isospin mode exchanged between the two-baryon states 31 and 24, as in fig. 
acts like an effective interaction. In this appendix, we derive an expression for 
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the effective spin-isospin interaction M (34, 12), wliicli consists of excfiange of all the 
spin-isospin modes present in our RPA approximation. This effective interaction then 
appears naturally in the calculations of the A width and A-cross sections, as presented 
in sec. |]|. M(34, 12) is of the form 



M''*(34, 12) = #'*(31) [?9''*(24)]* M''*(34, 12) 



Lti 



fit/ 



(132) 



where the spin-isospin matrix elements, ^?''*(31), are defined in eqs. (|60D and (|6T|). 

As before we will present two different expressions for M(34, 12). The first orig- 
inates from the summation of non-interacting Green's functions according to the 
diagrams in fig. |I], while the second expression originates from the RPA expansion in 
eq. (|57|). The first expression becomes 



M'(34, 12) 
M*(34, 12) 



D^Fyqlsi34, 12) + F^g'^^i^A, 12) 



fn fTT 
J 31 J 24 

;D,FW^{3A, 12) + JM^F^d'eA^A, 12) . 



With 



hu = 63- ei = 62 - 64, q = P3-Pl=P2-Pi ^ 

the dressed pion, D^,, and p-meson, Dp, propagators are written 
with the polarization functions 



n.(g) 



F2 
\2 ^-K 



det 



{RffxN + {RtfxA + F'^ 



h Xn Xa 



(133) 
(134) 

(135) 

(136) 

(137) 



and 



np(g) 



-icqi? 



F' ml 



^ j^N fNN 

/det V fNN, 



XN + Ri 



A fNA 



fN 



Xa 



NA, 



+ Fg Xn Xa] 



(138) 



The form factors F.,^ p and Fg, depend on the transferred energy and momentum q = 
{huj,cq), and are defined in eqs. (^) and (|T3p. The relativistic corrections appearing 
in d), (I), (0), and |12D are here denoted i?", with 



a 



(139) 



The index i denotes that center-of-mass energy y/s originates from an internal vertex, 
and is thus approximately taken 



Si{q) = {rriNC^ + hivf - {cqf 



(140) 
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Furthermore is the Ntt center-of-mass momentum at an internal vertex, and is 
approximated by 



rriNC 



2^* -9- (141) 
rrLNC^ + nuj 

The susceptibihties Xa are defined from the Lindhard functions in eqs. ( |121D and 



Xa(g) = M, 



In 



Na 



a 



Mr 



4 a = N 
16/9 a = A 



In the expressions of -Dvr.p there also appears the renormalization factor 

/det(g) = 1 + ^3 9'nn Xn + Fg g'^^ XA + Fg 6g[ xn Xa , 
and we have used the short hand notation 

\2 



^93 



J 



9nn9aa ~ {9'na) 
(Ra) 9'nn + (Rn) 9'aa - '^RnRa9'na 



(Ra) 9nn 



NA 



(Rn) 9'aa 



NA, 



'f_ 

Jn 



NN 



'^R^nR\9na 



NA 



NN , 



J NNJ N 



NA 



The quantity g^g consists of four terms 

g2g(34, 12) = to(34, 12) + toi(34, 12) + tio(34, 12) + tn(34, 12) 



with 



(142) 



(143) 



(144) 
(145) 

.(146) 



(147) 



to(34,12) 
tio(34, 12) 

toi(34,12) 
til (34, 12) 



R R QSI- 924 5 



dct 



9'n,24 Rf Xn + 9 AM Rf Xa-^ 



+ bg'^RfFlxNXA 



tio(21,43) 



5^ 

n 



dot 



148) 



(149) 
(150) 



9n,z\ 9n,2a K^i ) Xn + 9a.3i 9a,2a K^i ) Xa+ 



+ (5'iV,31 9'a,2A + fi'A,31 fi'7V,24) Rf Rf XA XN+ 



+ ^91 {9'n,3i Ri + 5'7v,24 Ri )Ri Fg Xn Xa 
+ 59'i {9'a,31 RT + ^?k,24 Rf)Rt F; xn Xl 
+ {5g[rRrRrF^x'Nxl\ , (151) 

where the relativistic correction at the external n-jk or p-jk vertex is taken as 



2m,7cC^ 



(152) 
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with 

sA^) = el- icp,r ■ (153) 

In this paper we only take into account vr or p vertices with either two nucleons or one 
nucleon and one A, and there is no ambiguity of which mass to substitute for rrijk. 
The notation means that the relativistic correction originates from an internal 
vertex, but the baryon mass appearing in Rf^ is determined from the baryons at the 
external jk vertex, 

inw = ^"^ff . (154) 



The A^vr {N p) center-of-mass momentum in the external vr-jfc (p-jk) vertex is given 

by 



Finally, the quantity (7gg(34, 12) is given by 

^eff(34, 12) = 24 - ^[9'n,u9'n,2aXN + ^k,31^A,24XA + ^31,24 ^d'l XTVXa] • (156) 

Jdet 

Alternatively, M(34, 12) can be expressed using an expansion in RPA eigenstates, 
M''*(34,12) = [ ^'"(3^''^)^';^^''^) h^'\?>l-,v)h^'\2A-v) \ 

2X2^9^^,12 ■ (157) 



The factors h^'^{jk; v) are motivated in eq. ( p5| ) and can, using the RPA solution in 
appendix ^ be explicitly expressed as 



a=N,A 



where the quantities Xa and are defined in appendix ^ and 



a,jk/ \ I J N,aJ N,jk t-i2/ „\ /1 c;n\ 

V (^, Q) = 9a,jk ^^^^2)2 ^9 (^' 9) • (159) 
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mjv = 940 MeV/c^ 
rriA = 1230 MeV/c^ 

= 140 MeV/c2 
rrip = 770 MeV/c^ 


9'nn = 0.9 
g'NA = 0.38 
^kA = 0.35 
Ag = 1.5 GeV 


fNN = 1-0 

Ina = 2.2 

/Ia = 

A^ = 1.0 GeV 


Inn — 6.2 

fNA = 10.5 
/^A = 

Ap = 1.5 GeV 


po = 0.153 fm-^ 


Va-Vn^ 25.0 MeV 


m^ = m;v/[l + 0.4049(p/po)] 



Table 1: Parameter values used in the numerical calculations. 
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Figure 1: 

Diagrammatic illustration of the microscopic structure of the Green's function G^^^ 
for a spin-isospin mode. A given spin-isospin mode (wiggly line) consists of dressed 
IT OT p mesons (long-dashed) and dressed baryon-hole excitations. The center row 
shows how a dressed meson is built from the corresponding free meson (short dashed) 
and dressed baryon-hole excitations. The double (solid-dashed) line in the dressed 
baryon-hole bubble can represent either a nucleon or a A isobar, as illustrated in the 
bottom row, where a single solid line represents a nucleon and a double solid line 
represents a A isobar. 

Figure 2: 

Diagrammatic representation of the effective spin-isospin interaction M(34, 12). 



Figure 3: 

Diagrammatic representations of the A self energy Sa (left-hand side) and the partial 
width (right-hand side). 



Figure 4: 

Differential cross sections in vacuum for the process p+p A+++n, in the pp center- 
of-mass system. In (a) is shown da/dcos{9) for y/s = 2.314 GeV, while (6) shows 
da/dt for i/i = 2.513 GeV. The sohd curve is the full cross section, da / d cos^Ocm) , the 
dashed curves are the direct and exchange contributions from the spin-longitudinal 
channel, and the dash-dotted curves are the contributions from the spin-transverse 
channel. There are also contributions from mixed direct and exchange terms; these 
are included in the full cross section but not displayed since they are small. The data 

I, but have here been estimated from figs. 5 and 8 



points originate from refs. ^1 
of ref. 137 



Figure 5: 

The total cross section in vacuum for the process p-|-p — A+++n, in the pp center- 
of-mass system as a function of y/s. The solid curve is our calculation, the dashed is 
the parameterization from VerWest-Arndt |2^, and the dash-dotted line is a simple 
parameterization often used in BUU calculations |^ . The data points were estimated 
from fig. 2 of ref. 



23 and can be found in references therein. 



Figure 6: 

The quantity [F*^* - F|^^ + (5F^'^^'']/2, calculated at the density 0.75p^ and for two 
different values of the parameter g'^^; the A energy and momentum are related to 
the pion kinetic energy by eq. (^). The displayed quantity may be compared with 
the spreading potential in ref. p5|. The empirical points were originally determined 
in ref. [^, but have here been estimated from fig. 12 of ref. |T^. 
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Figure 7: 

The dispersion relations for the spin-isospin modes in infinite nuclear matter at normal 
nuclear density and zero temperature. Parts (a) and {b) show the real part of ^a;^ 
in the spin-longitudinal and spin-transverse channel, respectively. The corresponding 
imaginary parts are presented in (c) and (d). The non-collective modes are shown by 
solid curves, while collective modes are represented by either a dot-dashed curve (vfi 
or pi), a dot-dot-dashed curve or P2), or a dot-dot-dot-dashed curve (see text). As 
a reference, the free pion dispersion relation htUj^^q) = [(m^c^)^ -|- (cg)^]^/^ is included 
as a dotted curve. 



Figure 8: 

Squared amplitudes for the spin-isospin modes displayed in fig. 0a, for q = 300 MeV/c: 
the individual NN~^ components (a), the individual AA^~^ components (6), and the 
pion component together with the sum of all NN"^ and AA^~^ components (c). 
Moreover, for the lower collective mode vfi is shown the q dependence of the pion 
component and the total NN"'^ and AA^~^ components {d). 



Figure 9: 

Same as fig. |^, but for the density Pn = '^Pn temperature T = MeV in 

(a) and for p^v = Pn T = 25 MeV in (b). 



Figure 10: 

The total A width F^* calculated either with Faat-i = in eq. ( |122|) (a) or with Faat-i 
included self-consistently (6), for a variety of scenarios: piv = 0, T = (dotted), 
Pn = Pat, T = (solid), = 2p5^, T = (short-dashed), p^ = p%, T = 25 MeV 
(long-dashed), and Pn = '^P%, T = 25 MeV (dash-dotted). 



Figure 11: 

The total A width F^* and its partial contributions from different spin-isospin modes. 
The solid curve represents the total width, the long-dashed line is the contribution 
from the non-collective NN~^ modes, the short-dashed line is the contribution from 
the non- collective AA^^^ modes, the dot-dashed line is the contribution from the 
lower pionic mode vfi, and the dot-dot-dashed line is the contribution from the upper 
pionic mode tt2. The calculations have been made with Faat-i = in eq. ( |122|) and 
for: p^ = p^, T = (a), pn = 2p%, T = (6), pn = p%, T = 25 MeV (c), and 
Pat = 2p^, T = 25 MeV (d). The error bars indicate the estimated uncertainty 
associated with the classification procedure (see text). 



Figure 12: 

Similar to fig. |ll|, but with Faat-i included self-consistently in eq. (|2 
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Figure 13: 

Differential cross section da /dcos{d cm) for tlie process p+p^ + n in tlie nuclear 
medium, in tlie pp center-of-mass system at the energy ^/s = 2314 MeV, calculated 
either with F^tv-i = (a) or with F^^v-i = T^* (6), for the following scenarios: 
Pn = 0, T = (dotted), pn = p%, T = (solid), p^ = 2p%, T = (short-dashed). 



Pn 



T = 25 MeV (long-dashed), and pn = 2p%, T = 25 MeV (dash-dotted). 



Figure 14: 

Total cross section cr{\^) for the process p+p + n for either Faat-i = in 

eq. ( |122| ) (a) or with T^^^-i included self-consistently (6). The notation is the same 
as in fig. |T^. 



Figure 15: 

Same as in fig. but for the reverse process n -t-A 
m = 1230 MeV/c^, with either Faat-i = in eq. 
self-consistently (6). 



p + p, for a A with mass 
(a) or with Faat-i included 



Figure 16: 

Total cross section o^{^/s = m) for the process p+TTj 
self-consistently and for the following scenarios: pn - 
T = (solid), Pn = '^P%, T = (short-dashed), pn 



P% 



A++ with Faat-i included 
0, T = (dotted), PN = P%, 
T = 25 MeV (long-dashed). 



and Pn = '^p'n, T = 25 MeV (dash-dotted). Parts (a) and (6) are for the lower pionic 
mode TTi, while parts (c) and (d) are for the upper pionic mode 7r2. The open squares 
represent the empirical total cross section for vr+p scattering and are estimated from 
fig. 2.2 in ref. 0. The factor p^irn^) in eq. ( pBf ) was in (a) and (c) calculated from 
the free A width, while the total width in nuclear matter was used in (6) and {d). 



Figure 17: 



Same as in fig. |TT|a and |TT|c, but the intermediate nucleon in the A decay is not Pauli 
blocked. 
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